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Preface 


Working  Group  No.  10  of  the  AGARD  Fluid  Dynamics  Panel  (FDP)  was  formed  to  investigate  the  limits  of  boundary  layer 
methods,  both  the  integral  and  field  type  formulation,  for  the  calculations  of  three-dimensional  turbulent  separate  flows. 

Since  the  first  meeting  in  late  1985  various  difficulties  arose  with  regard  to  technical  and  organizational  problems  and  the 
present  report  reflects  this  in  some  respect.  However,  it  is  hoped  that  the  results  presented  give  the  reader  sufficient  inspiration 
to  appreciate  the  work  performed  by  the  various  members  of  the  Working  Group  and  their  co-workers. 

The  participants'  names  are  given  below. 


Le  groupe  de  travail  No.lO  du  Panel  AGARD  de  la  dynamique  des  fluides  a  ete  constitue  pour  examiner  les  limitations  des 
methodes  de  calcul  appliquees  a  la  couche  limite,  tant  a  formulation  du  type  integral  que  du  type  champ,  pour  le  ealcul  des 
ecoulements  turbulents  decodes  tridimensionnels. 

Depuis  la  premiere  reunion,  organisee  vers  la  fin  de  I'annee  1985,  un  certain  nombre  de  problemes  d’ordre  technique  et 
logistique  se  sont  poses  et  le  present  rapport  en  temoigne  dans  une  certaine  mesure.  Cependant,  il  est  a  esperer  que  les  resultats 
qui  y  sont  prwentes  permettront  aux  lecteurs  de  se  rendre  compte  de  I'amplitude  des  travaux  entrepris  par  les  differents 
membres  du  groupe  et  leurs  collaborateurs. 
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CHAPTER  1.  BACKGROUND  AND  OVERVIEW 


H.Yoshihan,  Office  of  Naval  Research,  Liaison  Office. 
Far  East,  APO  San  Francisco  96503-0007 


1.  INTRODUCTION 

Flow  separation  plays  a  dominant  role  In  the  high  lift  performance  of  combat 
aircraft.  Invariably  limiting  takeoff  and  transonic  maneuver  performance.  Such 
separations  arise  In  two  forms.  The  first  Is  In  the  form  of  free  shear  layer 
separations  which  for  example  arises  In  the  well  known  fashion  along  highly  swept, 
sharp-nosed  wings  at  large  angles  of  attack.  The  second  Is  the  bubble-type  separation 
where  the  separated  flow  Is  confined  within  the  separation  bubble  as  In  the  case  of 
shock-induced  separations  on  an  unswept  wing  In  the  transonic  regime.  To  treat  such 
viscous  vortical  flows.  It  Is  essential  that  properly  modeled  viscous  flow  equations 
be  used,  which  for  relevent  cases  will  be  turbulent. 

With  the  recent  availability  of  vector  supercomputers  with  their  large  speeds 
and  memories,  there  has  been  Impressive  progress  In  computations  using  the 
Reynolds-averaged  Navler/Stokes  equations.  Despite  use  of  sophisticated  algorithms 
programmed  In  effective  vector  formats,  the  computation  costs  have  been  so  large  as  to 
preclude  their  use  by  the  aerospace  community.  A  viable  affordable  alternative  Is 
needed  for  design  applications,  and  this  Is  provided  at  large  Reynolds  numbers  by  the 
boundary  layer  procedure,  either  In  the  differential  (field)  or  Integral  formulations. 

In  the  boundary  layer  method  the  problem  Is  divided  Into  two  parts,  the  boundary 
layer/wake  and  the  equivalent  Inviscid  flow.  In  the  direct  formulation  the  equivalent 
Inviscid  flow  problem  Is  first  solved  assuming  the  viscous  effects  to  be  negligible. 
The  resultant  surface  pressure  and  flow  direction  distributions  are  then  Inputted  into 
the  turbulent  equations  simplified  by  the  boundary  layer  approximation.  Calculation 
of  the  boundary  layer  equations  then  yields  the  familiar  displacement  thickness  or  Its 
equivalent  the  surface  transpiration  velocity.  Either  of  the  latter  Is  then  Inputted 
Into  the  equivalent  Inviscid  flow  to  update  the  configuration  tangency  condition  to 
reflect  the  viscous  displacement  effects.  The  above  cycle  Is  then  repeated  until 
convergence  of  the  Iterative  cycle  Is  achieved.  In  the  absence  of  separation,  the 
above  direct  formulation  was  found  to  be  convergent  with  the  solution  In  most  cases 
comparing  closely  with  experiments. 

In  the  direct  mode  described  above,  the  field  equations  are  of  mixed 
parabolic/hyperbolic  type,  whereas  the  Integral  equations  In  most  formulations  are 
fully  hyperbolic.  In  the  presence  of  flow  separation,  the  above  Iterative  procedure 
In  the  direct  mode  ceased  to  converge  when  a  marching  direction,  successful  In  the 
direct  mode,  was  used  to  solve  the  boundary  layer  problem.  The  primary  cause  of  the 
difficulty  was  an  Instability  of  the  boundary  layer  procedure  (In  both  the  field  and 
Integral  methods)  due  to  the  finite  difference  marching  not  being  In  the  proper 
downwind  direction.  Moreover,  In  the  field  method  a  ’saddle-point”  type 
111 -conditioning  has  been  experienced  for  near-unswept  separations. 

Additionally,  In  the  Integral  formulation,  when  the  form  factor  H  ("delta 
star/theta")  assumed  values  In  the  range  2-3,  a  singularity  arose  due  to  the  form 
factor  modeling,  halting  the  marching.  Since  the  form  factor  Is  In  most  cases  less 
than  the  ’singular  value’  for  free  shear  laye,  separations,  the  above  difficulty  was 
confined  to  the  more  severe  bubble  type  separations,  for  example,  arising  In  the 
transonic  case.  It  Is  not  clear  whether  a  form  factor  relation,  free  of  the 
singularity,  can  be  evolved  using  Improved  velocity  profiles.  For  this  purpose, 
appropriate  Navler/Stokes  solutions  should  be  examined  for  guidance. 

To  circumvent  the  problems  described  above  for  separated  flows,  the  boundary 
layer  problem,  both  In  the  field  and  Integral  forms,  was  posed  In  one  of  several 
Inverse  forms  where  two  of  the  boundary  layer  variables  were  used  as  Inputs  Instead  of 
the  surface  pressure  and  flow  direction.  In  cases  where  the  Inverse  formulation 
permitted  a  marching  (cases  not  partially  elliptic),  a  reorientation  of  the  ’downwind’ 
direction  then  permitted  the  use  of  the  convenient  marching  used  In  the  direct  mode. 
Additionally,  the  form  factor  singularity  In  the  Integra)  formulation  was  removed.  A 
difficulty  that  then  arose  with  the  Inverse  formulation  was  the  loss  of  the 
Input/output  compatibility  between  the  boundary  layer  and  equivalent  Inviscid  flow 
problems  which  then  precluded  the  classical  Iterative  procedure  of  the  direct  mode. 
In  the  Inverse  mode,  the  output  from  both  problems  was  the  surface  pressure  and  flow 
direction  which  In  general  would  not  be  the  same.  An  update  procedure  was  then  required 
to  determine  Improved  Inputs  to  both  problems  as  a  function  of  the  pressure/flow 
direction  mismatch.  Simplified  local  update  procedures  successful  In  the  planar 
problem  (see  for  example  LeBalleur,  Ref.  1),  were  found  to  be  Inadequate  In  the  3D 
cases  when  significant  separation  was  present  (Ref.  2). 

The  simplifications  due  to  the  boundary  layer  approximation  are  1. 
simplification  and  stiffness  reduction  of  the  equations,  and  2.  the  decoupling  of  the 
Inviscid  and  viscid  portions  of  the  flow.  The  same  simplifications  can  be  achieved 
without  Incurring  the  above  update  problem  by  patching  the  boundary  layer  flow  to  the 
outer  Inviscid  flow  using  a  multiblock  mesh.  Here  the  multiblock  mesh  from  the 
NASA/Ames  transonic  N/S  (TNS)  code  (Ref.  3)  might  be  used  In  which  the  Inner  blocks 


eabed  the  conffguratfon  with  its  boundary  layer  and  wake  and  the  outer  blocks  contain 
the  Inviscid  portion  of  the  flow.  In  the  Inner  blocks  the  unsteady  boundary  layer 
nethod  developed  by  Van  Dal sem  and  Steger  (Ref.  4)  would  be  used.  An  overlapping  of 
the  nesh  blocks  then  provides  the  franework  to  Impose  flow  continuity.  Another 

sinllar  approach  of  pronise  Is  the  Fortified  N/S  method  proposed  by  Van  Dalsem  and 
Steger  (Ref.  5).  It  avoids  the  difficulties  of  the  boundary  layer  approach  by 

overlapping  the  boundary  layer  and  the  outer  flow  and  effecting  the  coupling  through 
the  novel  use  of  a  source  term. 

In  summary,  the  situation  for  computing  separated  boundary  layer  flows  Is  as 

follows;  In  the  case  of  the  field  method,  one  must  march  In  the  proper  downwind 
direction  at  each  level  of  the  boundary  layer.  For  large  turning  of  the  flow  within 
the  boundary  layer  as  would  usually  arise  In  separated  flows,  care  must  be  exercised 
In  the  choice  of  the  proper  finite  difference  marching  direction.  To  ease  this  task. 
Van  Dalsem  and  Steger  (Ref.  4)  used  the  unsteady  formulation,  properly  biasing  the 

transverse  spatial  differences.  For  the  case  of  the  Integral  method,  there  Is  not 
only  a  need  to  march  In  the  proper  downwind  direction,  but  the  form  factor  singularity 
must  be  dispensed  with,  either  by  posing  the  problem  In  the  Inverse  mode,  or  In  a  more 
difficult  fashion,  by  employing  a  compatibility  condition  to  neutralire  the 
singularity  (Ref.  2).  In  the  inverse  mode,  the  difficulty  then  arises  of  a  need  of  an 
yet-to-be  evolved  update  procedure  to  permit  Iterations  between  the  boundary  layer  and 
equiva'ient  Inviscid  flows.  The  use  of  a  flow  patching  with  a  multiblock  mesh  and  the 
use  of  the  Fortified  H/S  method  Is  then  suggested  as  means  to  avoid  this  difficulty. 

2.  FLUID  DYNAMICS  PANEL  WORKING  GROUP  10 

In  view  of  the  prospects  of  treating  significant  separated  flows  with  greatly 
reduced  computing  costs  relative  to  the  Navler/Stokes  method,  the  AGARD  Fluid  Dynamics 
Panel  (FDP)  established  Working  Group  (WG)  10  “3D  Viscous  Flows-Boundary  Layer  Limit" 
with  the  following  Objective  and  Scope. 

Dbjecti ve 

The  primary  objective  Is  to  develop  procedures  to  calculate  compressible 
turbulent  separated  flows  over  lifting  bodies  and  wings  In  the  high  Reynolds  Number 
limit  to  be  carried  out  In  two  steps: 

1.  Develop  techniques  to  solve  the  boundary  layer  equations  for  separated  flows 
prescribing  appropriate  direct  or  Inverse  Input  functions. 

2.  Develop  techniques  to  couple  the  resulting  boundary  layer  problem  with  the 
equivalent  Inviscid  flow  problem. 

Scope 


The  WG  will  be  confined  to  steady  turbulent  flows  using  both  the  field  and 
Integral  boundary  layer  equations.  In  the  former,  the  emphasis  will  be  on  algebraic 
turbulence  models  though  differential  equation  models  may  also  be  considered.  The 
primary  goal  here  will  not  be  on  the  validation  of  the  turbulence  model,  but  on  the 
development  of  the  mechanics  of  the  boundary  layer  procedure  and  Its  coupling  with  the 
Inviscid  flow  method. 


Two  test  cases  to  be  considered  are  the  DFVLR  6:1  prolate  spheroid  and  the 
Lockheed/NASA-Ames  Wing  C  Fighter  Wing  for  which  experimental  data  are  available.  The 
Navler/Stokes  solutions  will  also  be  used  as  a  comparison  data  base. 


There  was  great  Interest  among  the  NATO  member  countries  In  the  WG  10,  but  there 
was  Initial  concern  as  to  how  to  approach  the  tasks  with  the  wide  disparity  In  the 
stage  of  development  In  the  various  countries.  It  was  therefore  agreed  upon  that  the 
WG  would  function  as  a  clinic  whereby  the  participants  would  exercise  their  methods  at 
their  respective  stages  and  share  their  experiences,  particularly  with  respect  to  the 
numerical  difficulties  encountered  and  their  cures.  This  cooperative  spirit  would  be 
carried  out,  respecting  the  proprietary  requisites  of  the  participating  groups. 
Members  of  the  WG  10  were  as  follows: 


Mr.  J.  C.  LeBalleur  Dr.  H.  U.  Meier 

DNERA-Chatlllon  DFVLR-Gottingen 


Dr.  P.  D.  Smith 
RAE-Farnborough 


Dr.  B.  van  den  Berg 
NLR-Amsterdam 


Prof.  H.  Norstrud  Dr.  J.  Steger 

University  of  Trondheim  NASA-Ames 


Prof.  T.  Cebeci 
CSU-Long  Beach 


Prof.  M.  Onorato  Dr.  H.  W.  Stock 

Polltecnico  d1  Torino  Dornler  GmbH 


Dr.  C.  Gleyzes 
ONERA-Cert 


Mr.  A.  Pagano  Or.  H.  Voshihara 

British  Aero$pace-F11 ton  Boeing  Company 


Prof.  E.  H.  Hirschel  Dr.  A  Panaras 

MBB-Ottobrun  Athens  (NASA'Ames) 


Mr.  D.  Humphreys  Prof.  L.  Roberts 

FFA-Stockholm  Stanford  University 
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The  first  chairaan  of  the  WG  was  Dr.  H.  Yoshihara,  who  in  turn  was  succeeded  by  Dr.  A. 
Panaras  and  Prof.  H.  Norstrud.  (This  change  in  the  chairmanship  was  necessitated  by 
the  AGARD  policy  requiring  the  chairaan  to  be  a  aeaber  of  the  Fluid  Dynamics  Panel.) 
The  term  of  the  WG  was  from  September  1985  to  June  1988.  WG  meetings  were  held  over 
two  days  during  the  week  preceding  the  FOP  business  meeting.  Meeting  were  held  as 
follows:  Fall  1985  (AGARD  HQ,  Paris):  Spring  1986  (Polltecnico  d1  Torino,  Torino); 
Fall  1986  (Stanford  University,  Palo  Alto):  Spring  1987  (Landon);  Fall  1987  (Naples); 
and  Spring  1988  (Lisbon).  Tasks  were  carried  out  by  the  following  Committees; 


1 . 

2. 

3. 

4. 

5. 

The  remainder 


Data  Base:  Chairman,  Dr.  B.  van  den  Berg 

a.  6:1  Prolate  Spheroid,  Or.  H.  U.  Meier 

b.  Wing  C,  Dr.  A.  Panaras 

c.  Navier/Stokes  Sojtlons,  Or.  J.  Steger 

Integral  Methods:  Chairman,  Dr.  P.  D.  Smith 
Field  Methods;  Chairman,  Or.  T.  Cebeci 

Inviscid/Viscid  Interaction:  Chairman,  Mr.  J.  C.  Le  Balleur 

Coordination  with  Eurovisc  Committee  and  "Stanford  II" 
Committee  (Prof.  G.  Lllley):  Mr.  D.  Humphreys 

of  the  report  is  organized  as  follows: 


Chapter  2.  Physical  Aspects  of  30  Separated  Flows 
Or.  B.  van  den  Berg 


Chapter  3.  Theoretical  Formulation 

Mr.  J.  C.  LeBalleur*  and  Prof.  Dr.  E.  H.  Hirschel 


Chapter  4.  Numerical  Procedures 

a.  Integral  Methods 

Dr.  P.  D.  Smith*.  Or.  C.  Gleyzes,  and  Dr.  H.  Stock 

b.  Field  Methods 

Mr.  0.  Humphreys*.  Dr.  T.  Cebeci,  and  Dr.  fl.  van  den  Berg 

Chapter  5.  Data  Base  for  Test  Cases 

a.  Experimental  Cases 

Or.  H.  U.  Meier*,  and  Dr.  H.  Yoshihara 

b.  Navier/Stokes  Calculations 

Or.  0.  Steger*,  and  Or.  A.  Panaras 

Chapter  6.  Summary  of  Results  and  Computational  Guidelines 
Or.  L.  Roberts*,  and  Mr.  D.  Humphreys 

a.  Integral  Methods 

Dr.  P.  Smith*,  Mr.  J.  C.  LeBalleur,  Dr.  C.  Gleyzes,  and 
Or.  H.  Stock 

b.  Field  Methods 

Dr.  T.  Cebeci*,  Dr.  B.  van  den  Berg,  Prof.  Dr.  E.  H. 

Hirschel,  Mr.  D.  Humphreys,  and  Or.  F.  Menter 

Chapter  7.  Conclusions. 

Prof.  Or.  H.  Norstrud*,  and  Or.  H.  Yoshihara 

The  asterisks  denote  the  principal  authors. 
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CHAPTER  2 

PHYSICAL  ASPECTS  OP  3D  SEPARATED  BOUNDARY  LAYER  FLOWS 

B.vui  den  Berg,  National  Aerospace  Laboratory.  NLR, 

Anthony  Fokkerweg  2, 1059  CM  Amsterdam,  The  Netherlands 

In  this  chapter  tvo  8ub;)ectB  vlll  be  dlacuaaed  sore  specifically.  In  the  first  place  a  description 
will  be  given  of  30  separated  flow  patterna  in  general.  This  Beans  that  attention  will  be  focussed  on  the 
topological  features  of  the  flow.  It  is  important  to  understand  these  features  for  the  calculation  process 
to  be  followed  and  for  the  interpretation  of  the  results.  As  practical  flows  are  generally  turbulent  and  an 
exact  description  of  turbulence  is  far  beyond  reach  at  thla  moment,  an  approximate  mathematical  model  for 
the  turbulence  terms  in  the  equations  has  to  be  employed  In  calculations.  These  empirical  models  are  at 
least  partly  based  on  general  physical  notions  and  they  are  the  second  subject  to  be  discussed  in  this 
chapter. 

When  referring  to  calculations,  field  methods  for  the  viscous  flow  have  been  assumed  here,  as  they 
relate  more  simply  to  the  physics.  Integral  methods,  which  solve  the  boundary  layer  equations  In  a  global 
way,  actually  do  a  comparable  job,  but  the  turbulence  tsodel  is  hidden  in  global  empirical  assumptions.  Also 
Che  discussion  on  the  topology  of  3D  separated  flows  is  not  basically  different  for  integral  methods, 
including  the  properties  near  separation  (Coustelx  &  Houdevtlle  1981). 

2.L  Topological  features  of  separated  flows 

It  toay  be  useful  to  start  recapitulating  shortly  the  features  of  separation  in  2D  stationary  flows. 
Separation  is  then  simply  identified  with  a  line  with  skin  friction  t  •  0,  Indicating  the  beginning  of  a 
reverse  flow  region,  which  la  Inaccessible  for  fluid  from  upstream.  Further  features  are  that  the  boundary 
layer  thickness  grows  generally  rapidly  at  and  downstream  of  the  separation  line  and  that  the  interaction 
between  the  viscous  and  invlscld  flow  becomes  large.  It  is  well  known  that  solutions  of  the.Aamlnar 
boundary  layer  equations  with  pressure  dlscribution  prescribed  produce  a  singularity,  t  ~  x^'  ,  near 
separation. Similar  behaviour  may  be  expected  for  turbulent  boundary  layer  calculations,  dependent  on  the 
turbulence  model.  However,  regular  solutions  can  be  obtained  when  prescribing  other  boundary  conditions 
such  as  the  boundary  layer  displacement  thickness  variation  (e,g.  Catherall  &  Mangier  1966).  These,  as  well 
as  full  Mavler-Stokes  solutions,  indicate  that  the  physical  flow  behaves  regularly,  t  "  x,  near 
separation. 

In  3D  flows  the  skin  friction  “  0  only  at  singular  points  and  not  generally  along  a  line  on  the 
surface.  Close  to  such  a  point  the  skin  friction  components  t  and  t  In  the  orthogonal  directions  x  and 
y  may  be  written  in  case  of  regular  behaviour:  *  .x  +  ^  .y  an^T  -a  .x  +  a  .y,  retaining  only 

the  loweat-order  terms  in  Che  series  expansion,  wtch  means  tha?  the  disctfsslo^f’Vlll  b^^restricted  here  to 

first-order  singularities.  Dependent  on  the  sign  and  magnitude  of  the  coefficients  a  .  a  ,  a  ,  a  , 

different  types  of  singular  points  are  obtained  <e.g.  Lighthlll  1963).  Figure  2.1  ahoi^  pattern^of  ftle 
skin  friction  lines  or  wall  streamlines  near  such  points  for  some  typical  cases.  If  J  «  a  .a 
a  >  0  an  inftnita  number  of  wall  streamlinea  pass  through  the  singular  point.  The  singular  poli^  Is 

called  a  node  in  this  case.  If  A  •  a  +  a  >  0  it  is  a  node  of  attachment  as  sketched  in  figure  2.1. a. 
When  changing  the  direction  of  the  arrowA^of  the  wall  streamlines  In  the  sketch  a  node  of  separa¬ 
tion,  A  <  0,  is  obtained.  Clearly  on  an  isolated  three-dimensional  body  at  least  one  node  of  attachment  and 
one  node  of  separation  must  be  present.  If  6/4  -  J  <  0  the  wall  streamlines  approach  the  singular  point 

asymptotically  in  spirals.  Such  a  node  is  called  a  focus,  aee  figure  2.1.6.  The  presence  of  a  focus  of 

»«paraclon  generally  merka  ch-  '<^velopment  of  a  vortex  in  the  flow.  If  J  <  0  the  wall  streamlines,  except 
two,  do  not  pass  through  the  singular  point.  Such  a  point  la  calxeu  a  Raddle  end  If  A  <  0  it  is  a  saddle 
point  of  separation  as  sketched  in  figure  2.l.c.  Finally  it  seems  useful  to  mention  in  this  short  survey  a 
well-known  topological  rule,  which  states  that  the  number  of  nodal  points  (nodes  and  foci)  exceeds  the 
number  of  saddle  points  by  two  on  bodies  with  no  multiply  connected  surfaces  (Lighthlll  1963). 

Nodes  are  either  sources  or  sinks  of  wall  streamlines,  as  contrasted  with  saddles.  The  special 

feature  of  a  saddle  is  that  it  partitions  the  wall  streamlines  coming  from  both  sides  towards  the  saddle. 
The  "wo  wall  streamlines  leaving  the  saddle  act,  at  least  locally,  as  a  barrier  between  wall  streamlines 
coming  from  opposite  directions.  An  exceptionally  efficient  barrier  occurs  in  the  degenerate  case  of  a  line 
with  saddle  points  of  separation.  Such  a  line  corresponds  to  what  is  called  a  separation  line  in  2D  flows. 
Analogously  one  might  call  the  wall  streamlines  leaving  a  saddle  point  separation  lin-’R  fu  3D  flow®.  More 
precisely  Lighthlll  (1963)  has  proposed  to  define  as  separation  lines  wall  streamlines,  which  Issue  from  a 

saddle  point  of  separation  and  disappear  in  a  node  of  separation.  It  should  be  noted,  however,  that  many  of 

the  phyalcal  features  attending  2D  separation,  as  sunarlzed  earlier,  do  not  apply  here.  The  differences 
will  be  discussed  in  some  detail  In  the  next  paragraphs. 

In  the  first  place  in  3D  flows  the  akin  friction  is  not  zero  except  at  the  saddle  point  itself. 
Further  it  is  not  generally  true  that  the  wall  streamlines  Issuing  from  a  saddle  point  act  as  a  complete 
barrier  between  flow  from  "upstream"  and  "downstream",  as  a  2D  separation  line  does.  This  is  Illustrated  in 
figure  2.2,  which  shows  some  simple  wall  flow  patterns  near  one  or  two  nodes  and  saddles.  It  Is  evident 
from  figures  2.2.b  and  c  chat  Che  wall  streamlines  issuing  from  the  saddle  points  and  disappearing  in  the 
foci,  do  not  constitute  a  complete  barrier.  Actually  the  wall  streamlines  directed  towards  the  saddle  point 
from  downstream  all  originate  upstream,  which  means  that  the  whole  surface  is  accessible  from  upstream. 
Figure  2. 2. a,  shoving  the  flow  pattern  for  a  saddle  with  a  node  downstream,  provides  an  example  where  the 
wall  streamlines  Issuing  from  ths  saddle  Indeed  constitute  a  complete  barrier.  Even  In  this  case,  however, 
this  holds  only  for  the  wall  streamlines  and  not  necessatily  for  any  other  streamline  away  from  the  wall, 
since  the  external  streamline  attaching  to  the  downstream  node  normally  does  not  originate  at  the  upstream 
saddle,  see  figure  2. 2. a.  A  closed  region  with  separated  fluid,  which  is  inaccessible  for  fluid  from 
upstream,  as  In  2D  separations,  does  not  exist  in  the  case  shown.  In  general  it  may  be  stated  that  closed 
separation  regions  probably  are  earlier  an  exception  than  the  rule  In  3D  flows  (Hunt  et  al.  1978).  The  feet 
that  3D  flow  separations,  due  to  the  extra  dloMnelon,  generally  tend  to  be  open  separations,  seems  not  to 
be  always  sufficiently  well  appreciated. 


s 

In  2D  flows*  separation  is  accompanied  with  rapid  boundary  layer  growth.  Such  a  rapid  growth  does  not 
necessarily  occur  in  3D  flows  near  wall  streamlines  issuing  from  saddle  points  of  separation.  To  evaluate 
the  growth,  the  continuity  equation  for  the  flow  near  the  surface  will  be  considered.  Assuming  U  "  t  .z» 
one  obtains  t  .z  .An  •  constant,  where  z  is  the  wall  distance  and  An  is  the  distance  between  two 
neighbouring  wall  streamlines.  In  2D  flows  An  *  constant,  but  t  ■»  0  at  separation,  so  that  the  wall 
distance  z  must  grow  very  rapidly  to  satisfy  continuity.  In  3D  flows  the  skin  friction  does  not  go  to  zero, 
except  at  the  singular  points,  and  rapid  boundary  layer  growth  along  special  wall  streamlines  is  not 
evident.  Rapid  growth  must  come  in  that  case  from  a  fast  decrease  of  An,  l.e.  from  a  strong  flow 
convergence.  Strong  flow  convergence  does  not  always  occur  along  wall  streamlines  issuing  from  saddle 
points.  On  the  other  hand  experimental  evidence  indicates  that  strong  flow  convergence  may  occur  along  wall 
streamlines,  which  do  not  issue  from  a  saddle  point,  and  «rhlch  do  not  distinguish  themselves  in  any 
fundamental  aspect  from  other  wall  streamlines.  Further  downstream  along  such  a  wail  streamline  a  very 
thick  boundary  layer  will  have  developed  with  many  of  the  physical  features  that  are  generally  associated 
with  separation.  This  has  led  Wang  (1976)  to  call  this  an  "open  separation".  The  terminology  is  confusing, 
however,  since  most  3D  separations  are  open,  as  argued  In  the  preceding  paragraph.  The  important  conclusion 
is  that  separation-like  flow  features  may  start  to  develop  along  a  wall  streamline  at  a  position,  which  is 
not  well  defined  and  Is  not  a  singular  point,  while  along  wall  streamlines  issuing  from  a  saddle  point 
separation-like  features  may  not  become  very  evident  at  all. 

Since  the  boundary  layer  thickness  increases  rapidly  at  and  downstream  of  a  2D  separation  line,  the 
viscous  effect  on  the  Inviscld  flow  will  be  large.  On  the  other  hand  small  changes  of  the  pressure  gradient 
imposed  by  the  inviscld  flow  on  the  viscous  flow  may  have  large  effects  as  the  velocities  in  2D  separation 
regions  are  small.  This  is  a  physical  explanation  for  the  strong  vlscous-invlscld  Interaction  in  these 
conditions.  Because  rapid  boundary  layer  growth  does  not  always  occur  in  3D  flows  along  Che  lines  which  may 
be  defined  as  separation  lines,  and  because  the  velocities  are  not  always  small  in  3D  separation  regions, 
strong  viscous-inviscid  Interaction  does  not  necessarily  occur  here  on  physical  groundb.  A  mathematical 
reason  for  strong  vlscous-lnviscid  interaction  may  be  singular  behaviour  of  the  solution  near  the 
separation  line,  like  in  2D  flows  when  Che  boundary  layer  equations  are  solved  with  the  pressure 
distribution  prescribed.  The  behaviour  of  the  solution  of  Che  3D  boundary  layer  equations  near  special 
lines  like  wall  streamlines  Issuing  from  saddle  points,  has  not  been  investigated  thus  far  in  detail  for 
the  general  case.  For  quasi-two-dimensional  flows  some  conclusions  can  be  drawn,  however,  from  a  simple 
analysis  to  be  discussed  hereafter. 

In  the  analysis  laminar  boundary  layers  will  be  assumed,  but  the  results  are  likely  to  hold  also  for 
turbulent  flows  as  with  Che  normal  turbulence  models  the  equations  are  similar  to  a  fairly  high  order  of 
approximation  close  Co  the  wall,  which  is  Che  region  of  interest  here.  For  quasi-two-dimensional  laminar 
boundary  layers,  as  occur  on  infinite  swept  wings,  the  so-called  independence  principle  holds  (e.g. 
Schllchtlng  1979).  The  Independence  principle  states  that  the  velocity  component  normal  to  the  lines  of 
invarlancy  (e.g.  Che  infinite  swept  wing  leading  edge)  can  be  solved  independently  from  the  velocity 
component  parallel  to  these  lines.  This  means  chat  chat  part  of  the  solution  is  identical  to  the  2D 
solution.  It  can  be  concluded,  therefore  chat  in  Chat  case  wict^^e  pressure  distribution  prescribed  the 
skin  friction  component  normal  to  Che  separation  line  t  ~  x  near  separation,  as  in  2D.  The  other 
component  t  •  constant,  ^cnaequently  the  skin  friction^  wall  streamlines  touch  In  cusp-llke  manner  the 
separation  ^ne,  which  is  then  an  envelope  of  wall  streamlines.  The  solution  thus  obtained  is  singular  at 
the  separation  line.  Like  in  2D  flows,  however,  regular  solutions  may  be  obtained  when  employing  more 
appropriate  solution  techniques,  applying  ocher  boundary  conditions  or  solving  Che  full  Navler-Stokes 
equations.  In  such  calculations  the  skin  friction  near  separation  will  vary  as  t  ^  x,  t  •  constant, 
which  means  chat  the  separation  line  is  an  asymptote  of  wall  streamlines. 

The  question,  whether  a  3D  separation  line.  Issuing  from  a  saddle  point  and  disappearing  In  a  node. 
Is  an  envelope  or  an  asymptote  of  wall  streamlines,  has  received  continual  attention  In  the  literature 
(e.g.  Maskell  1955,  Llghthill  1963,  Brown  &  Stewartson  1969,  Wang  1976,  Hornung  &  Perry  1964).  If  the 
foregoing  qua8i-cwo-di.&tei.£ional  analysis  may  be  generalized,  it  is  clear  Chat  a  separation  line  is  not  an 
envelope  of  wall  streamlines,  neither  in  physlcel  reality  nor  in  proper  calculations.  This  conclusion 
should  not  be  surprising,  since  singularities,  as  occur  at  ruspa,  can  not  be  expected  to  exist  In  reality. 
The  occurrence  of  cusps  at  separation  in  boundary  layer  calculations  Is  due  to  Improper  boundary  conditions 
and  may  be  dependent  also  on  the  marching  direction.  In  fully  3D  flows  the  wall  streamlines  probably  do  not 
really  approach  Che  separation  line  asymptotically,  but  all  touch  Che  separation  line  at  the  downstream 
node,  where  it  terminates.  Surface  flow  visualizations  see-*  *‘o  suggest  Chat  separation  lines  are  envelopes 
oi  wall  streamlines,  but  obviously  rigorous  conclusions  can  never  be  drawn  from  such  obseivetioub. 

It  will  be  evident  from  the  foregoing  chat  the  definition  of  the  separetlon  line  in  3D  flows  poses 
problems.  An  unambiguous  definition  is  that  of  Llghthill:  a  separation  line  is  a  wall  streamline  Issuing 
from  a  saddle  point  of  separation  and  disappearing  in  a  node  of  sepsratlon.  Many  of  the  physical  flow 
features  associated  with  2D  separations  may  not  occur  along  such  vail  streamlines,  however.  Also,  as  these 
lines  are  no  envelopes  of  wall  streamlines,  they  may  not  be  diatlngulahable  from  any  other  wall  streamline 
apart  from  their  special  origin.  At  Che  same  time  vail  streamlines,  which  do  not  originate  at  a  saddle 
point,  may  show  many  of  the  physical  features  generally  associated  with  separation.  Since  in  practice  one 
normally  wants  to  Include  the  latter  case  in  the  separated  flows,  several  other  definitions  for  a  3D 
separation  line  have  been  considered,  none  of  them  satisfactory  (see  e.g.  Vang  1976).  Definitions  based  on 
the  breakdown  of  the  boundary  layer  assumptions  or  rapid  boundary  layer  growth  are  inadequate  as  these  do 
no  provide  precise  criteria.  Whether  singular  solutions  of  the  boundary  layer  equations  occur  along  a 
separation  line  depends  on  the  solution  technique.  Insccessiblllty  for  upstream  fluid  also  does  not  lead  to 
a  satisfactory  discrimination  as  mostly  ths  vholt  flow  is  accessible  in  the  3D  case.  It  seems  that  the 
separation  line  concept  adopted  from  2D  flows  is  not  simply  transfsrsble  to  3D  flows. 

If  the  term  separation  is  yet  employed  for  3D  flows,  as  is  usual  practica,  it  is  profitable  to 
distinguish  two  Cypts  of  ssparation  llnss:  wall  streamlines  issuing  from  saddls  points  and  lines  of 
convergence  not  involving  singular  points.  Thera  is  considsrsbls  diversity  in  the  terminology  for  the  two 
types  of  ssparation  in  ths  lltsraturs,  s.g.  bubble  vs  frse-shesr-laycr  separation  (Maskell  1955),  closed  vs 
open  ssparation  (Wang  1976),  global  vs  local  ssparation  (Tobak  &  Peake  1982),  while  also  ths  term  cross¬ 
flow  ssparation  is  used  for  the  latter  type  (Chapman  1986).  Unfortunately  all  the  terms  proposed  are 
misleading  to  soma  extent  due  to  preconceived  ideas  about  the  flow  physics,  which  do  not  always  apply.  In 
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face  flov  Aft  further  downstreaa  behind  e.g.  a  saddle  and  a  node  (figure  2. 2. a)  say  be 
Indlatingulahabie  from  those  In  tvo  regions  of  flow  convergence  Involving  no  singular  points.  Therefore,  it 
seeu  mott  rational  Co  refer  to  both  types  of  separation  on  the  bisls  of  their  eleaentary  distinguishing 
feature.  This  suggest  the  use  of  teras  like  saddle  point  separation  and  gradual  separation  respectively. 
Since  in  the  latter  case  no  singular  points  are  involved  and  the  separation  characteristics  develop 
grsjually,  the  Cera  separation  Is  Inherently  not  precise  here. 

In  this  report  two  teat  cases  for  boundary  layer  calculations  are  considered.  These  include  the  3D 
turbulent  boundary  layer  flow  on  an  elllpaold  at  tO*  and  at  30*  angle  of  attack  and  a  transonic  swept  wing 
flow  with  a  local  separation.  The  flow  aeparatlona  on  the  elllpaold  at  angle  of  attack  are  typical  gradual 
separations  Involving  no  alngular  points.  Though  far  downatreaa  on  the  ellipsoid  saddle  point  separations 
■ay  occur,  cha  3D  separation  in  Che  region  of  interest  here,  l.e.  Che  region  where  boundary  layer 
calculations  have  been  perforaed,  la  of  Che  gradual  type  at  both  angles  of  attack  (Meier  eC  al  1983).  On 
Che  avape  wing,  ho%rcver,  the  surface  flow  visualizations  (Keener  1985)  indicate  the  presence  of  a  saddle 
point  separation.  Actually  a  saddle  and  a  focuS  are  cXeuly  recognizable  on  the  wing  upper  surface  at  the 
test  case  conditions  (M^  ■  0.85,  a  ■  5*.  Re  ■  6.8  x  10”).  The  separated  wall  streamline  pattern  Is  basi¬ 
cally  similar  to  that  sketched  in  figure  2.2^b. 


2.2  Turbulence  aodelllng  of  separated  flows 

As  a  universal  model  for  Che  turbulent  shear  stresses  in  the  Reynolds-averaged  flow  equations  does 
not  yet  exist,  and  probably  will  never  do,  seml-emplrlcal  models  have  to  be  applied  In  calculation  methods. 
These  seml-emplrlcal  models  have  generally  a  restricted  range  of  applicability.  This  leads  co  turbulence 
models  valid  for  specific  flow  zones.  In  this  section  turbulence  models  for  separated  flows  will  be 
reviewed  briefly.  Actually  the  emphasis  will  be  on  the  complexities  and  the  consequent  unsatisfactory  state 
of  affairs  In  turbulence  modelling,  notably  for  separated  flows. 

It  may  be  useful  to  start  again,  as  In  section  2.1,  with  a  discussion  of  the  situation  for  2D 
separated  flows.  When  the  boundary  layer  is  not  separated,  the  flow  properties  In  the  wall  region  determine 
to  a  large  extent  the  2D  turbulent  boundary  layer  development.  Since  the  flow  In  the  wall  region  seems  to 
behave  reasonably  well  universally  (as  appears  from  the  wide  validity  of  the  law  of  the  wall),  simple 
turbulence  models  describing  at  least  the  wall  region  correctly  produce  often  satisfactory  results.  The 
altuaClon  becomes  different,  however,  when  the  flow  separates  from  the  wall.  The  boundary  layer  becomes 
then  a  free  shear  layer  developing  above  a  recirculation  region  In  which  the  velocities  are  comparatively 
small.  Aa  the  stabilizing  effect  of  Che  wall  is  absent,  the  turbulence  level  In  the  free  shear  layer  tends 
to  Increase  substantially.  When  the  separated  flow  reattaches,  the  free  shear  layer  is  splltted  Into  two 
parts,  one  part  flowing  backwards  Into  the  recirculation  region.  This  will  affect  especially  the 
development  of  the  large  turbulent  eddies.  Clearly  substantial  changes  In  turbulence  level  and  structure 
occur  in  case  of  a  2D  separation  (and  reatcachment) .  As  these  changes  cake  place  generally  comparatively 
fast,  turbulence  history  effects  may  be  expected  to  play  an  important  role.  In  addition,  2D  separated  flows 
are  often  unsteady  and  Che  low  frequency  fluctuations,  t^lch  probably  do  not  depend  on  local  conditions 
only,  may  affect  the  turbulence  structure  and  can  produce  a  non-negllglble  direct  contribution  to  the 
turbulent  ahear  scressea  (e.g.  Bradahaw  1978).  It  will  be  clear  that  complex  turbulence  models  are  needed 
to  describe  the  phenomena  in  such  flows  more  or  less  satisfactory.  Actually  the  accuracy  of  the  turbulence 
models  available  up  to  now  for  2D  separation  regions  Is  Inadequate  (e.g.  Kline  et  al  1963,  Delery  &  Marvin 
1986). 


As  argued  In  section  2.1,  the  physical  features  ettendlng  a  3D  separation  do  not  necessarily  resemble 
very  closely  those  of  s  20  separation.  For  Instance,  If  there  Is  a  recirculation  region  downstream  of  a  3D 
separation  line,  the  velocities  need  not  to  be  small  In  Chat  region  as  there  may  be  a  significant  velocity 
component  parallel  co  the  separation  line.  Moreover,  often  3D  separations  do  not  feature  a  distinct 
recirculation  region,  the  separation  region  being  open.  It  Is  also  not  possible  generally  to  distinguish 
downstream  of  a  3D  separation  line  clearly  a  shear  layer  sway  from  the  wall,  as  in  2D  separations.  Actually 
cha  flow  near  and  downstream  of  a  3D  separation  line  may  not  be  very  different  from  a  normal  attached  3D 
boundary  layer  flov.  It  will  be  evident  from  Che  foregoing  chat  the  substantial  changes  in  turbulence  level 
end  structure  attending  2D  separations,  do  not  necessarily  occur  In  the  same  way  in  3D  separations.  On  the 
ocher  hand  the  flov  la  generally  strongly  three-dimensional  near  3D  separations  and  the  usual  problems 
associated  with  turbulence  modelling  for  3D  boundary  layers  apply  here  consequently. 

It  la  generally  accepted  now  Chat  simple  scalar  eddy  viscosity  turbulence  models  are  Inadequate  for 
3D  boundary  layers  as  a  considerable  eddy  viscosity  non-lsotropy  appears  co  exist.  Some  experimental  values 
of  ths  ratio  of  the  crosswise  to  the  strcamwlse  eddy  viscosity  N  In  three-dimensional  boundary  layers  are 
collected  In  figure  2.3.  The  plot  Includes  Che  recent  experimental  data  from  Anderson  and  Eaton  (1967), 
which  again  show  chat  cha  ratio  may  attain  values  substantially  less  chan  one.  Though  the  non-isotropy 
of  the  eddy  vlscoalty  has  attracted  most  attention,  the  comparatively  low  magnitude  of  the  eddy  viscosity 
in  3D  boundary  laysrs  according  co  expsrlments  is  at  least  ss  Important  to  take  into  account  for  accurate 
turbulence  modelling.  Figure  2.6  shows  Che  considerable  decrease  of  the  eddy  viscosity  magnitude  measured 
In  s  boundary  layer  developing  from  a  2D  condition  Co  a  3D  separation  (Elsenaar  &  Boelsma  1976). 

In  general  txperlmenca  Indicate  that  the  turbulence  level  decreases  when  the  flov  becomes  three- 
dimensional.  Varloua  explanatlona  have  put  forward.  It  has  been  proposed  that  due  to  the  flov  Chree- 
dlmensionallty  the  large  turbulent  eddies  are  hindered  In  their  development  as  they  will  topple  and  that 
consequently  turbulence  activity  decreases  (Bradshaw  6  Pontlkos  1985).  Further,  history  effects  may  be 
expected  to  be  reaponslblc  for  at  least  part  of  Che  effects  found  (Van  den  Berg  1962).  Three-dimensional 
boundary  layers  Involve  generally  also  flov  convergence  or  divergence,  which  are  known  to  lead  to  a  sub¬ 
stantial  turbulence  level  dsersase  or  increase  respectively  (s.g.  Pstel  6  Bsek  1987).  Due  to  the  extra 
dimension,  additional  phenomena  which  affect  the  turbulence  properties  occur  In  3D  flows  and  complicate 
accurate  turbulence  modelling. 

It  may  be  concluded  that  turbulence  modelling  in  3D  aeparatlon  regions  poses  different,  chough  not 
necessarily  more  difficult  problems  than  in  2D  separation  regions.  In  both  cases  ths  problems  arc  tar  from 
solved,  however,  and  it  looks  like  that  in  the  not  too  distant  future  the  probably  still  existing 
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shorcconlngs  of  che  turbulence  models  will  be  the  principal  factor  limiting  the  accuracy  of  boundary  layer 
calculations,  notably  in  separation  regions. 


References 

Anderson,  S.D. ,  Eaton, J.K.,  1987,  "An  experimental  investigation  of  pressure  driven  three-dimen::ional 

turbulent  boundary  layers".  Stanford  University,  Therm.  Div.  Rep.  MD-49. 

Berg,  B.  van  den,  1982,  "Some  notes  on  three-dimensional  turbulent  boundary  layer  data  and  turbulence 
modelling' . 

lUTAM  Symposium,  Berlin,  Eds:  H.H.  Fernholz,  E.  Krause. 

Bradshaw,  P. ,  1978,  "Structure  of  turbulence  in  complex  flows".  AGARD  LS  No.  94. 

Bradshaw,  P. ,  Pontlkos,  N.S.,  1985,  "Measurements  in  the  turbulent  boundary  layer  on  an  Infinite  swept 
wing". 

J.  Fluid  Mech.  159.  105. 

Brown,  S.N.,  Stewartson,  K.,  1969,  "Laminar  separation".  Annual  Review  of  Fluid  Mech..  Vol.  1,  45. 

Catherall.  D. ,  Mangier,  K.W. ,  1966,  "The  integration  of  the  two-dimensional  laminar  boundary  layer 
equations  past  che  point  of  vanishing  skin  friction”  J.  Fluid  Mech.  26,  163.  1966. 

Chapman,  G.T. ,  1986,  "Topological  classification  of  flow  separation  on  three-dimensional  bodies". 

AIAA  Paper  86-0485. 

Coustelx,  J.,  Houdeville,  R.,  1981,  "Singularities  in  three-dimensional  turbulent  boundary  layer 
calculations  and  separation  phenomena".  AlAA  J.  19,  8,  976. 

Dechow,  R.  ,  1977,  "Hlttlere  Ceschwlndigkelt  und  Reynoldsscher  Spannungstensor  in  der  driedlmensionalen 

turbulenten  Wandgrenzschicht  vor  elnea  stehenden  Zylinder". 

Dlss.  Un.  Karlsruhe,  Strdmungsmech.  &  Stromungsmachinen  Heft  21,  1977. 

Delery,  J.,  Marvin,  J.G. ,  1986,  "Shock-wave  boundary  layer  interactions".  ACARDograph  No.  280. 

EXsenaar,  A.,  Boelsma,  S.H.,  1974,  "Measurements  of  the  Reynolds  stress  censor  in  a  three-dimensional 
turbulent  boundary  layer  under  infinite  swept  wing  conditions".  NLR  TR  74096  U. 

Hornung,  H, ,  Perry,  A.E.,  1984,  "Some  aspects  of  three-dimensional  separation".  Part  I:  "Stream  surface 
bifurcations".  Zeitschrift  Flugwise.  &  Welcraumforschung  8,2,77.  1984. 

Hunt.  J.C.R.,  Abell,  C.J.,  Peterka,  J.A.,  Woo,  K. ,  1978,  "Ktnemacical  studies  of  the  flows  around  free  or 
surface-mounted  ooaCsclesi  applying  topology  to  flov  visualization'* ,  J.  Fluid  Mech.  86,  179. 

Keener,  E.R.,  1985,  "Pressure  distribution  measurements  on  a  transonic  low-aspect  ratio  wing". 

NASA  TM  86603. 

Kline,  A.J.,  Cantwell,  B.J.,  Lllley,  C.M.,  1982,  "Complex  turbulent  flows". 

The  1980-81  Stanford  Conf . ,  Vol.  2. 

Lighchill,  M.J,,  1963,  "Attachment  and  separation  in  three-dimensional  flow".  In:  L.  Rosenhead  (ed.): 
"Laminar  boundary  layers",  Oxford  Un.  Press. 

Maskell,  E.C.,  1955,  "Flov  separation  In  three  dimensions".  RaE  Aero  Report  2565. 

Meier,  H.U. ,  Kreplln,  K.P.,  Vollmers,  H.  1983,  "Development  of  boundary  layers  and  separation  patterns  on  a 
body  of  revolution  at  incidence",  Second  Symp.  "Numerical  and  Physical  Aspects  of  Aerodynamic  Flows", 

Ed. :  T.  Cebecl. 

Patel,  V.C.,  Baek,  J.H.,  1987,  "Boundary  layers  In  planes  of  symmetry;  Part  1:  Experiments  In  turbulent 
flow".  AIAA  J.  25,  4,  550. 

Pontlkos,  N.S.,  1982,  "The  structure  of  three-dimensional  turbulent  boundary  layers". 

Ph.D.  Thesis,  Imperial  College,  London. 

Schllchting,  H. ,  1979,  "Boundary  layer  theory".  McCraw  -  Hill. 

Tobak,  M.  ,  Peake,  D.J.,  1982,  "Topology  of  three-dimensional  separated  flows".' 

Annual  Review  of  Fluid  Mech.,  Vol.  14,61. 

Wang,  K.C.,  1976,  "Separation  of  three-dimensional  flow".  Martin  Marietta  Labs  TR-76-54  C. 


node  of 

ATTACHMENT 


NORMAL  VALUE 


U) 


CHAPTER  3.  THEORETICAL  FORMULATION 
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The  Navier-Stokes  equations  are  believed  to  describe  any  phenomenon  of  continuum 
flow,  even  turbulence.  Discrete  numerical  solutions  of  the  Navier-Stokes  equations  are 
available  today,  but,  except  for  some  special  cases,  actually  the  Reynolds-averaged 
equations  are  solved  for  non-leuninar  flow.  The  special  cases  concern  large-eddy  or 
direct  simulation  of  turbulent  flow,  which  is  possible  on  a  research  level  now  (see  for 
instance  Ref.  1). 

Although  the  computer  power  has  grown  dramatically  in  the  last  decade,  solutions 
of  the  Reynolds-averaged  Navier-Stokes  eqx'ations  are  much  too  expensive  to  be  applied 
in  regular  aerodynamic  design  work,  Ref.  2.  In  addition  there  are  still  many  not  suf¬ 
ficiently  solved  problems  with  discrete  numerical  methods,  and  the  big  question  marks 
which  have  to  be  put  to  our  abilities  to  model  transition  laminar-turbulent  and  tur¬ 
bulence  in  separating  shear  flows. 

The  aerodynamic  design  engineer  traditionally  employs  a  whole  palette  of  flow 
models  and  computation  methods  in  his  work.  Depending  on  the  given  design  problem  he 
distinguishes : 

-  inviscid  flow, 

-  inviscid  flow  with  weak  interacting  viscous  effects, 

-  inviscid  flow  with  strong  interacting  viscous  effects, 

-  fully  viscous  flow, 

which  usually  appear,  although  in  different  areas,  and  with  different  importance,  on 
every  aircraft. 

Another  way  to  classify  flow  fields  is  by  distinguishing 

-  attached  flow, 

-  separated  flow. 

There  is  a  problem  with  the  concept  of  separation,  because  as  it  usually  is  used 
it  stems  from  two-dimensional  flows,  where  it  simply  is  defined  by  vanishing  surface 
shear.  Actually  a  more  general  view  should  be  adopted.  Refs.  3  and  4,  Whether  the 
boundary  layer  leaves  a  sharp  edge  (flow-off  separation  at  the  trailing  edge  of  a  wing, 
or  at  the  leading  edge  of  a  delta  wing  at  larger  angle  of  attack) ,  or  gets  separated 
at  the  body  surface  (free-surface  or  squeeze-off  separation,  the  latter  especially  in 
three-dinnsnsional  flow  situations),  in  each  case  the  following  effects  appear: 

-  strong  interaction  of  the  separating  boundary  layer  with  the  inviscid  flow, 
which  leads  to  a  viscosity  induced  pressure  drag, 

-  convective  transport  of  vorticity  away  from  the  body  surface  in  wakes  (vortex 

sheets)  and  vortices  (the  induced  drag  of  a  finite-span  lifting  wing  is  re¬ 

lated  to  this  phenomenon) . 

If  one  accepts  this  more  general  view  it  becomes  evident  at  once  that  separated 
flow  fields  in  some  cases  can  be  treated  by  means  of  purely  inviscid  models,  like  po¬ 
tential  flow  models,  see  for  instance  Ref.  5,  or  by  Euler  models.  Ref.  4,  see  also  the 
very  important  Ref.  6.  The  typical  application  considered  in  these  references  is  that 
to  delta  wings  with  sharp  leading  edges  at  angles  of  attack,  where  leading-edge  separa¬ 
tion  appears  and  hence  lee-side  vortices. 

The  scope  of  the  present  study  is  limited  to  free-surface  or  squeeze-off  separa¬ 
tion,  where  the  location  of  separation  is  not  fixed  by  sharp  edges.  The  question  is 

whether  and  how  far  flow  fields  with  such  properties  can  be  handled  with  some  inviscid 

flow  -  boundary-layer  approach.  The  basic  principle  which  allows  such  an  approach  is 
the  locality  principle.  Ref.  3, 

Experimental  and  theoretical/numerical  results  from  flowfield  investigations  on 
many  configurations  suggest  that  a  change  in  body  shape  or  in  the  flow  configuration, 
for  instance  by  separation,  is  felt  only  locally  and  downstream  of  that  location.  This 
is  the  reason  why  in  general  a  boundary- layer  computation  yields  good  results  compared 
to  experimental  results,  although  usually  only  inviscid  solutions  without  modelling  of 
the  separation  phenomena  are  applied. 

The  fact  that  also  the  location  of  primary  separation  lines  usually  can  be  pre¬ 
dicted  to  a  good  degree  of  accuracy  is  connected  with  this  principle.  The  global  inter¬ 
action  due  to  the  separating  boundary  layer  obviously  is  weak  in  most  cases.  The  local 
interaction,  however,  laecause  of  the  elliptic  properties  of  the  flow  near  a  separation 
line  or  point,  finally  makes  it  impossible  to  predict  exactly  the  location  of  separation 
in  the  frame  of  boundary- layer  theory. 


The  locality  principle  must  be  used  with  care.  Of  course,  the  flow  is  changed 
upstream  too,  for  instance  because  of  the  elliptic  property  of  subsonic  flow  fields. 
These  changes  are  small  as  can  be  seen.  They  can  be  significant,  if  for  instance 
the  wake  of  the  body  carries  kinematically  active  vorticity.  Then  a  global  inter¬ 
action  can  occur  which  may  not  be  significant  for  the  boundary-layer  development, 
but  for  the  global  forces  on  the  body.  An  example  for  this  is  the  induced  drag  of 
wings . 


The  formulation  of  boundary-layer  computation  problems  on  arbitrary  confi¬ 
gurations  poses  no  dif f icultiea.  Ref.  7.  Of  course,  boundary-layer  theory  must  be 
valid.  It  is  interesting  to  note  that  obviously  first-order  theory  rather  early 
becomes  incorrect  if  the  flow  is  decelerated,  even  if  the  boundary-layer  thickness, 
and  the  longitudinal  surface  curvature  is  not  large.  Ref.  8.  Second-order  boundary- 
layer  methods,  see  for  instance  Ref.  9,  allow  to  treat  thick  boundary-layers  on  curved 
surfaces,  where  centrifugal  pressure  corrections  must  be  made,  and  where  the  off- 
surface  metric  roust  be  tcdcen  into  account.  Of  course,  then  the  displacement  proper¬ 
ties  of  the  boundary- layer  must  be  imposed  on  the  inviscid  flow  by  means  of  the 
equivalent  inviscid  source  distribution  (Ref.  7),  whichis  done  routinely  in  zonal 
solutions.  Ref.  8.  Even  hypersonic  boundary  layers  with  entropy-layer  swallowing 
can  be  treated  with  second-order  boundary- layer  theory.  Ref.  10. 

The  coupling  of  a  boundary- layer  solution  with  an  inviscid  solution  in  order  to  com¬ 
pute  the  flow  in  a  separation  region  or  to  obtain  the  global  interaction,  poses  consider¬ 
able  problems.  In  two-dimensional  flow  cases  inverse  or  interactive  approaches  are 
often  used.  In  general  it  is  possible  to  extend  these  to  three-dimensional  flow  cases 
only  if  the  features  of  the  separated  flow  allow  this.  If  the  flow  exhibits  a  comp¬ 
licated  separation  topology,  then  the  Reynolds-averaged  Navier-Stokes  equations  have 
to  be  employed  globally  or  at  least,  like  in  the  zonal  solution  of  Ref.  8,  locally. 

The  boundary- layer  equations  are  of  parabolic  type.  Ref.  3  (or  in  the  Integral 
formulation  of  hyperbolic  type),  and  pose  an  initial  value/boundary-value  problem.  The 
stre^unlines  of  the  boundary  layer  are  characteristics,  and  therefore  the  marching  dir¬ 
ection  of  the  boundary-layer  method,  and  hence  the  grid  must  be  oriented  more  or  less 
in  main  streeun  direction.  A  complete  orientation  along  the  external  inviscid  stream¬ 
lines  (streamline  coordinates)  in  general  proves  to  be  too  complicated  to  handle,  be¬ 
cause  the  external  inviscid  streamlines  converge  or  diverge  in  parts  of  a  configuration. 
Then  streamlines  have  to  be  deleted  or  newly  to  be  started  (see  for  instance  Ref.  11). 

A  complete  configuration  orientation  of  the  grid,  like  in  the  cross-section  oriented 
grids  for  fuselages  or  the  percent-line  oriented  grids  on  wings.  Refs.  7  and  12,  also 
can  have  drawbacks.  Nose  or  attachment  regions,  also  the  regions  ahead  of  primary 
separation  may  pose  problems  with  such  grids.  An  interesting  approach  to  combine 
streamline  oriented  and  configuration  oriented  grids  in  so-called  hybrid  grids  was 
made  in  Ref.  13.  There  selected  Inviscid  streamlines  form  a  skeleton  in  which  then  a 
kind  of  skeleton  cross-section  grid  is  constructed.  In  this  way  the  advantages  of  both 
types  are  combined,  while  the  disadvantages  are  mostly  avoided. 

The  problem  of  initial  data  is  closely  connected  with  the  grid-orientation  problem. 
Although  one-dimensional  Initial  data  can  be  constructed  in  stagnation  points  or  on 
attachment  lines  in  chord  sections,  Ref.  7,  in  the  latter  case  in  the  infinite-swept  wing 
mode  (or  in  the  locally  infinite-swept  wing  mode,  where  the  correct  metric  properties  of 
the  wing  are  taken  into  account),  the  procedures  arc  quite  cumbersome.  Usually  therefore 
approximate  initial  data  are  employed  somewhat  downstream  of  the  actual  stagnation  point 
or  of  the  actual  attachment  line.  This  is  possible,  because  the  initial  data  very  fast 
loose  their  influence  on  the  downstream  solution  in  favour  of  the  local  boundary  data 
if  the  flow  is  sufficiently  accelerated.  Ref.  7.  However,  if  approximate  initial  data 
are  put  in  regions  of  decelerated  flow,  large  errors  may  build  up  in  the  complete  down¬ 
stream  solution,  unless  there  again  the  flow  is  accelerated,  as  for  instance  the  flow 
past  automobiles  behind  the  motor  hood.  Ref.  14. 
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CHAPTER  4.  NUMERICAL  PROCEDURES 

4.1  Integral  methods 


PJ).Siiiith,  AES  Division,  Royal  Aerospace  Establishment, 
Famborough,  Hants  GU 14  6TD,  United  Kingdom 


4.1.1  The  integral  equations 


Integral  prediction  methods  for  three-dimensional  turbulent  boundary  layers  are  most  readily  des¬ 
cribed  in  a  streamline  coordinate  system  (s^n)  which  consists  of  two  families  of  mutually  orthogonal 
curves  on  the  body  surface.  One  family  is  formed  by  the  projections,  onto  the  surface,  of  streamlines 
just  external  to  the  boundary  layer.  The  direction  (s)  of  an  external  streamline  is  called  the  streamvise 
direction  and  boundary-layer  flow  in  a  direction  (n)  normal  to  an  external  streamline  and  parallel  to  the 
surface  is  called  crossflow.  For  general  applicability  the  equations  in  the  s,  n  streamline  coordinate 
system  are  transformed  to  a  general  non-orthogonal  curvilinear  coordinate  system;  x,  y  in  the  surface  z 
normal  to  the  surface  with  an  angle  X(x,y)  between  the  x  and  y  directions  and  an  angle  a(x,y) 
between  the  x  and  s  directions. 


The  integral  equations  in  streamline  coordinates  are 
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8  momentum  integral: 
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Entrainment  or  continuity  integral: 
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Here  1]  is  the  external  velocity  in  the  streamline  direction,  H  denotes  the  Hach  number  corresponding 


to 


fs 


and 


are  the  skin  friction  components  in  the  streamvise  and  crosswise  directions  res¬ 


pectively.  6  is  the  boundary-layer  thickness, 
and  the  integral  thicknesses  are  defined  as: 


Cp  the  entrainment  coefficient,  H 
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In  addition  we  define 
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The  terms  K  and  K  are  the  geodesic  curvatures  of  the  curves  s  *  constant  (normals)  and 
8  n 

n  ”  constant  (streamlines)  respectively.  General  expressions  for  and  in  terms  of  the  angle  a 

aiul  the  geometry  of  the  x,y  coordinate  system  may  be  found  in  Reference  1 .  ia  a  measure  of  the  rate 

of  divergence  of  the  external  streamlines  whilst  if  the  z  component  of  the  vorticity  of  the  external 
flow  ia  zero  ■  1/U^  (3U^/3n  )  and  some  simplification  of  equations  (1)  to  (3)  results. 

Here  general  versions  of  equations  (I)  to  (3)  which  make  allowance  for  pressure  gradients  normal  to 
the  surface  and  include  terms  resulting  from  the  Reynolds  normal  stresses  may  be  found  in  References  2, 

3  and  13. 


4.1.2  Particular  methods  of  solution 

The  atajor  differences  between  the  integral  calculation  methods  used  by  the  members  of  the  Working 
Group  occur  in  the  assumptions  made  regarding  the  entrainment  coefficient  and  the  forms  assumed  for 

Che  atreamwise  and  crossflow  velocity  profiles.  These  profiles  are  needed  in  order  to  derive  relation¬ 
ships  between  the  integral  thicknesses  (equation  (4))  thus  reducing  the  number  of  unknowns.  The  major 
assumptions  made  are  detailed  below. 


14 


4. 1 .2. 1  RAE  method^ 

The  method  is  essentially  that  of  Smith,  Reference  4,  but  since  extended  to  allow  for  history  effects 
upon  the  turbulence  structure  (Reference  5)  by  calculating  the  entrainment  coefficient  by  means  of  the  lag 
equation  of  Green,  Weeks  and  Brooman  (Reference  6)  applied  along  the  external  streamlines.  Power  law  velo¬ 
city  profiles  are  used  in  the  streaowise  direction  and  Mager  (Reference  7)  profiles  in  the  crossflow  direc¬ 
tion.  To  improve  the  prediction  of  the  streamwise  flow  the  Uj  *  Hj(H)  relationship  used  is  not  that 

which  would  be  given  by  the  power  law  velocity  profiles  but  instead  the  »  H|(H)  relationship  as  given 

by  Lock's  (Reference  8)  equation  (105)  is  used.  This  includes  allowances  for  the  effects  of  low  Reynolds 
number  and  of  strong  departures  from  equilibrium.  In  addition  it  is  necessary  to  assume  a  skin  friction 
law  for  and  that  of  Green  et  al  (Reference  6)  is  used.  The  dependent  variables  of  the  method  are 

0|j,  H  and  B=  tan“'  . 

4. 1.2. 2  RAE  method  using  Cross's^  profiles 

Here  the  simple  but  fairly  restricted  power  law  and  Mager  profiles  of  Smith's  method  are  replaced  by 
the  more  satisfactory  but  complicated  profiles  suggested  by  Cross  (Reference  9). 
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is  the  standard  law  of  the  wall 
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F(c.x)  -  sin^jn?) 


Is  a  modification  of  the  Coles  wake  function  (and  allows  for  strong  departures  from  equilibrium)  and  C 
and  D  are  chosen  to  give  U/U^  ■  1,  V  ■  0  at  the  edge  of  the  boundary  layer  (4*1).  So  that 

C  •  1  -  cos  e  (In(U^Re^)  ♦  2.131  . 

D  •  (C  -  1)  tan  8  .  (7) 


The  exponents  Xg  as  given  by  Lock's  (Reference  8)  equation  (132)  and  give  improved  rep¬ 

resentation  of  profiles  with  large  values  of  6  •  With  this  f^ily  of  profiles  all  the  integral  thickness 

ratios  5./4,  6../5  can  be  expressed  as  functions  of  U  ,  Re.  and  6  and  no  additional  H, (H)  rela- 

1  IJ  TO  i 

tionship  and  skin  friction  law  are  required.  The  dependent  variables  of  the  method  are  5,  and  6  . 

4 . 1 . 2 . 3  ONERA  method 

This  uses  the  velocity  profiles  suggested  by  Le  Balleur^**®  He  replaces  the  wake  functions 

and  of  equation  (5)  with  the  single  function 

F(4,C*)  - 


where  C*  is  a  quasi-linear  function  of  di/d  ,  based  upon  the  computation  of  the  2D  mixing  layer,  which 
is  non  zero  only  for  grossly  separated  flows  (B  >  90^).  The  dependent  variables  of  the  method  are  6  , 
5|/5  and  ^2^^  *  entrairanent  coefficient,  ,  is  given  by  an  analytic  expression  deduced  from  the 

velocity  profiles  with  a  mixiog  length  assumption,  and  is  then  based  upon  equilibrium  flows  properties. 

At  present,  two  transport  integral  equations  (Reference  12)  for  non-equilibrium  effects  can  be  used  only 
in  two-dimensional  flows.  Minor  corrections  are  added  to  at  low  Reynolds  number  and  low  shape  para¬ 

meter.  No  additional  H| (H)  relationship  and  skin  friction  lav  are  required. 


4 . 1 . 2 . 4  ONBRA/CBRT  method 


Here  following  Cousteix  (References  15  and  16)  the  relationships  between  the  integral  thicknesses 
are  obtained  from  an  analysis  of  similarity  solutions.  These  similarity  solutions  are  constructed  by 
using  a  mixing  length  model  which  has  been  modified  to  account  for  the  fact  that  the  turbulent  shear 
stress  is  not  aligned  with  the  mean  velocity  gradient  3q/9z  .  These  effects  are  represented  by  a 
factor  T  where  T  ■  tanCyx  -  Y)/t«n(Yg  “  y)»  Y*  Yg  and  yt  are  the  directions  of  the  velocity,  the 
velocity  gradient  and  the  turbulent  shear  stress  respectively.  T  is  a  constant  for  a  particular  calcula¬ 
tion.  The  dependent  variables  of  the  method  are  ^||t  ^21  ~  their  equivalents  in  the 

general  x,y  coordinate  system).  The  cncrsinmenC  coefficient  is  derived  from  the  similsrity  solutions 
and  once  again  no  additional  Hj (R)  relationship  or  skin  friction  law  are  required. 
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4.1.3  Nusoerlcal  methods 

All  the  integral  methods  of  section  4.2  use  the  general,  non-orthogonal  coordinate  system  (x,y)  on 
the  body  surface.  The  results  of  inserting  the  assumptions  detailed  in  section  4.2  into  equations  (t)  to 
(3)  and  then  transforming  to  the  x.y  system  may  be  written  schematically  in  matrix  form  as 

aF.  ap. 

A..  -rJ- ♦  B..  -  C.  ,  (i. 3  -  1,2,3)  (8) 

ij  Sx  ij  3y  1 

where  is  the  vector  of  unknowns  (eg  6,  H,  6  for  Smith's  method)  and  the  and  are 

functions  of  the  unknowns,  the  geometry  of  the  x.y  aystcm  and  the  angle  a  .  The  right-hand  sides 
are  in  addition  functions  of  the  external  velocity  which  is  assumed  to  be  knotm.  For  attached  flows 

the  system  (8)  is  hyperbolic  and  can  be  solved  as  an  initial  value  problem  by  a  marching  procedure.  The 
bounding  characteristics  of  these  equations  lie  roughly  in  the  directions  of  the  external  and  surface 
streamlines,  ie  inclined  at  angles  a  and  a  -*■  6  respectively  to  the  x  axis. 

The  numerical  methods  of  &Dith,  RAE.  and  Gleyzes,  ONERA/CERT,  are  similar  in  that  both  are  explicit 
and  approximate  the  y  derivatives  of  equation  (6)  by  finite  differences  (backward  if  a  and  a  ■*>  B  are 
of  the  same  sign,  central  if  a  and  a  +  6  are  of  opposite  sign)  and  then  integrate  the  resulting  sys¬ 
tems  of  ordinary  differential  equations  with  respect  to  x  by  standard  methods  (two  step  Euler,  Smith; 
fourth  order  Runge  Rutta,  Gleyzes).  In  both  methods  the  x  step  is  chosen  to  satisfy  both  accuracy  and 
stability  (CFL)  requirements. 

Le  Balleur  et  al,  uae  the  second  order  explicit  KacCormack  scheme,  involving  central  differences 
for  the  y  derivatives,  «^ich  they  baV6 found  to  b6  more  robust  than  the  above  method.  In  addition  once 
either  of  the  characteristic  angles,  a  or  a  6  ,  become  large  x  marching  may  become  impossible  and 
Le  Balleur  has  developed  a  "Multi-Zonal-Marching"  method  in  which,  where  necessary,  marching  in  bothax 
anda  y  directions  is  employed,  referonoe  10. 

Usually  it  is  necessary  to  supply  boundary  conditions  at  the  lateral  (y)  boundaries  of  the  computa¬ 
tion  if  any  fluid  is  entering  the  COTputatlonal  domain  through  these  boundaries.  In  the  case  of  the 
ellipsoid,  this  is  not  necessary  since  the  flow  may  either  be  considered  to  be  periodic  in  y  ,  or  the 
boundaries  may  be  assumed  to  be  planes  of  synmetry. 
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4.2  Field  methods 


D.a.  Humphreys,  B.  van  den  Berg  and  T.  Cebeci 


A  non^orthogonal  curvilinear  coordinate  eyetee  <x,y,a).  useful  for  dlecuaeion  pur- 
poses,  has  been  defined  in  section  4a. 1  above.  Before  describing  the  individual  prograas 
though,  because  one  of  then  is  based  on  rather  unconventional  assuaptions,  it  eight  be 
appropriate  to  recall  here  a  little  of  the  background  to  the  subject. 

The  systea  of  boundary  layer  differential  equations  solved  In  field  aethods  is  a 
subset  of  the  exact  Reynolds  equations.  The  aethod  of  deriving  the  turbulent  boundary 
layer  equations,  although  well  supported  by  experimental  evidence,  is  necessarily  heuri* 
Stic  when  the  turbulence  ie  represented  only  by  empirical  foraulae.  For  laminar  flow 
matheaatically  rigorous  procedures  can  be  applied  and  in  that  case  lead  to  the  same 
simplifications  as  given  by  the  heuristic  aethod. 

The  essence  of  the  boundary  layer  concept  is  probably  contained  in  the  observation 
that,  with  the  very  high  Reynolds  nuabers  fortunately  prevailing  in  aircraft  aerodyna¬ 
mics.  diffusion  parallel  to  the  wall  tangent  plans  is  quite  negligible  and  exactly  zero 
at  the  surface.  The  rigorous  arguments  indicate  that  after  suitable  normalising,  and 
supposing  the  Reynolds  number  to  be  large  and  the  wall  curvature  to  be  eaall,  the  terse 
which  are  less  than  order  unity  srs  dominated  by  terse  of  order  either  IC4  or  3*,  where  K 
stands  for  the  maximum  curvature  of  any  local  perpendicular  section  through  the  wall  and 
4  is  the  local  boundary  layer  thicknese.  The  first  order  Prandtl  equations  consist  of 
the  0(1)  terms  alone  and  the  eecond  order  equations  retain  the  0(K3)  terse  but  reject 
those  of  order  4*1  this  imposes  s  consistency  limit  on  the  magnitude  of  the  peraitted 
wall  curvature*.  Both  aeta  of  boundary  layer  equatione  are  parabolic  and  in  steady  flow 
can  be  solved  numerically  by  space  marching,  given  the  requieite  flow-defining  initial 
ai4u  boundary  data.  In  the  second  order  theory  the  pressure  variation  normal  to  the  wall 
which  ie  associated  with  the  centrifugal  force  arising  through  wall  curvatura  amergas  as 
part  of  the  boundary  layer  calculation,  but  in  general  this  would  not  include  the  pres¬ 
sure  field  caused  by  streamline  curvature  near  separation  from  a  flat  wall  nor  that 
induead  for  Inatanee  by  longitudinal  vorticity. 

A  thorough  analysis  of  boundary  layer  calculation  methods  as  such  has  recently  been 
carried  out  in  another  context*.  Xt  appears  that  modarn  methods  are  superficially  very 
similar  and  in  ordsr  to  dif fsrsntiata  elsarly  ens  from'anothsr  it  Is  an  advantags  to  bs 
able  to  break  them  down  into  a  serias  of  logical  steps  and  eompars  how  each  slsmantary 
task  is  treatad.  A  similar  procedura  was  adopted  in  a  separate  document^  prepared  as  a 
atandardlasd  information  sourca  for  the  present  ehaper.  This  document,  together  with  the 
original  rsfareness  should  bs  consultsd  if  a  mors  complsts  prsscriptlon  of  s  cslculstion 
msthod  and  its  gsnsral  capabilitias  ars  sought.  Only  sssentials  will  bs  repestsd  hsrs, 
concsntrstlng  on  what  was  actually  dons  to  obtain  ths  results  appearing  in  this  report. 

There  ars  five  field  methods  to  describe  ss  in  chapter  is,  the  methods  will  bs 

referred  to  by  ths  organisation  where  they  were  devised:  AVA.  Odttingen;  Dac,  Long 
Beech;  FFA,  Stockholm;  NBB,  Mttnchen,  end  NLR,  Ameterdem.  The  HB8  method  solves  ths  se¬ 
cond  ordsr  boundary  layer  equations  whereas  ths  other  four  treat  ths  usual  first  ordsr 
Prandtl  equations.  Two  test  configurations  have  been  treated:  ths  ellipsoid  at  IQ^  Inci¬ 
dence,  which  sines  the  Mach  number  Is  only  about  0.10  may  bs  regarded  as  an  incompress¬ 
ible  flow,  and  the  NASA  wing  "C.  All  five  nethode  were  applied  to  ths  ellipsoid  flow. 
Only  ths  OAC  method  wee  used  for  wing  C  and  then  ths  energy  equation  was  solved  in  addi¬ 
tion  to  the  eompreesible  Prendtl  equations.  Towards  the  rear  of  the  ellipsoid  ths  boun¬ 
dary  layer  thicknese  becomes  quits  apprscisbis  in  relation  to  the  radius  of  the  body 
eroea  section  and  ths  question  arises:  should  the  metric  coefficients  bs  sllowsd  to  very 
with  distance  from  the  wall  as  ie  usual  with  programs  written  for  axl-symmetric  flow 
over  bodies  of  revolution?  zn  the  NBB  method  the  verietion  was  included,  but  neglected 
in  the  others.  Ths  principal  originator-ref ersness  where  general  descriptions  of  the 
field  methods  can  bs  found  srs  AVA  t4],  DAC  10,  6),  FFA  I?],  HBB  [8,  9,  20)  end 
NLR  (11). 

The  next  eabseetien  concerns  ths  turbulence  models.  Then  follows  s  definition  of 
ths  prlrcipml  numerical  methods  used  for  ths  computed  results.  Initial  end  boundary 
conditions  are  reserved  for  the  chapter  discussing  results,  beesuee  they  constitute  ths 
flow  definition  for  coaputationsi  purpoese  and  say  little  about  the  method,  on  the  other 
hand  a  note  on  computational  effort,  intimately  connected  with  numeiica,  is  appended  to 
the  present  chapter. 
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4.2.1  Turbulence  models 


The  turbulence  eodele  ueed  were  elgebreic.  ell  of  mixing  length  or  eddy  vleeoeity 
type.  Tor  ordinary  etteched  threedimeneionel  turbulent  boundary  layers  it  is  believed 
that  the  sophisticated  models  which  are  eometimee  proposed  for  practical  calculation 
methods*,  while  certainly  Implying  gross  Increases  in  computer  time,  cannot  in  general 


be  relied  upon  to  improve  the  quality  of  results  obtained.  However,  it  is  equally  car- 
tain  that  simple  extensions  to  three  dimensions  of  the  algebraic  models  developed  from 
obseryations  In  undisturbed  fully  developed  twodimensional  flow  are  not  generally  ade> 
quate^*.  The  flows  under  consideration  here  are  at  high  enough  Reynolds  numbers  for  well 
developed  turbulence  to  be  expected  but  some  of  the  other  factors  which  could  generate 
first  order  modifications  to  the  turbulence  are 


*20":  longitudinal  pressure  gradient 

longitudinal  wall  curvature 
transverse  wall  curvature 

*30":  streamline  convergence 

swirl 

eddy  viseoaity  anisotropy 

Of  these  the  influence  of  the  pressure  gradient  was  included  in  the  AVA,  DAC  and  NLR 
methods,  not  in  the  other  two.  None  of  the  remaining  effects  was  modelled.  The  equations 
used  may  be  given  as  follows,  assuming  streamline  coordinates. 


AVA  and 

DAC  methods  (Refs.  13,  S  and  21 

*  <  = 

“I  ■  ‘I'Kli)'* 

(4b:ll 

*e  <  *  = 

“T  ■ 

(4b:2) 

where 

-  axp 

(4b!2) 

Xj  -  0.4 

(4bt4) 

z'*'  ■  tifX/V 

(4b!5) 

»,  -  (1-11. e 

(4b:6) 

Ai*  -  26.0 

(4b!7) 

•  «■«/<> 

(4b;8) 

+  wu,  au, 

'1  u^«  ax 

(4b:9) 

•j  -  0.0168 

(4b;10) 

1-1  -  11. 0  +  8.8(z/a)*)*l 

(4b:ll) 

*1*  ■  I/O  <1-  5^1  <»*| 

(4b:12) 

■  u'  +  v* 

{4b!l3) 

Cm  is  the  total  wall  shear  stress  and  v  is  the  kinematic  viscosity  evaluated 
at  the  wall. 
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y>A  — thod  {iiAfs.  14  and  7] 
au^ 

S  <  X^t  l>,  -  t,«  5^ 

*C  <  *'  ■'r  ■ 

Hhara 

>2  "  Xiifl  -  xxp  -x*/K{*] 

•2  >  2/eO 

«2*  -  /s  f‘-  S;I  <** 


(4b;14) 

(4b:lS) 

(4b;16) 

(4b:17) 

(4b:18) 


Has  —thod  (Raf.  18] 

X  <  Xg!  I»,  -  •2*« 


<  x: 


T  “  «1>’2^3*“1 


whara 

rj  -  (1.0  +  9.a(0.3z/z..,)‘l'* 

83  “  l.**Bax 

Uj  «  M  aln(l,(u,/W)»/4] 


(4b:19) 

(4b:20) 


(4b:21) 

(4b;22) 

(4b:23) 


a  la  tha  aagnltuda  of  tha  vortlclty,  W  la  the  maxlaua  valua  of  the  quantity 
(eonaldarad  aa  a  function  at  z)  and  z.^^  Is  the  location  of  Its  naxlaum. 


WLR  nathod  (Rafa.  11  and  2] 


All  z; 


(4b!24) 


whara 


*3 

H2 

"2 

*2* 

+ 

•'2 

*1 

T(S) 

"1 


X2»ri  - 

axp  -z'*'N2/A2'''JTj 

(4b:25) 

0.41 

(4b;26} 

(1  -  10 

(4b:27) 

24. a 

(4b;28} 

MU.  au. 
**w 

(4b;29) 

T(»i) 

(4b:30} 

e"»  tanh  n 

(4b:31) 

X2Z 

0.11  4, 

(4b:32) 

x^  la  takan  In  tha  dlractlon  of  tha  wall  ahaar  atraaa  vactor  and  la  tha  valua  of  z 
*m«r«  ■  0.95  U0e 
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4.2.2  Numerical  methods 

The  meymptotie  dlecretleetlon  errors  of  the  celeuletion  sethode  have  the  orders  of  mag¬ 
nitude  shown  in  Table  1. 

Table  1  Orders  of  accuracy 


Method 

X  y  z 

AVA 

2  2  2 

DAG 

2  2  2 

FFA 

ill 

MSB 

2  2  1 

MLR 

112 

As  is  well  known,  the  orders  of  accuracy  give  a  measure  of  how  fast  the  local  discreti¬ 
sation  error  diminishes  as  the  step  length  goes  to  zero  and  do  not  in  themselves  carry 
any  Information  on  the  size  of  the  local  error  attached  to  any  individual  step  length 
taken  in  isolation.  First  order  methods  can  therefore  be  more  accurate  in  practice  than 
second  order  onea.  It  is  usually  an  unfortunate  fact  that  increasing  the  discretisation 
error  order  also  reduces  stability.  The  difficult  question  of  numerical  stability  will 
not  be  treated,  but  some  relevant  discussion  of  it  can  be  found  in  References  2  and  16, 
for  example. 

As  with  the  orders  of  accuracy,  the  integration  algorithms  are  conveniently  summa¬ 
rised  in  tabular  fora. 

Table  2  Inteoration  schemes 


Method 

parpandlcular 

tmngentiel 

AVA 

KB 

ZZ 

BAG 

KB 

GB  or  SB 

FFA 

IE 

LW 

NBB 

? 

? 

MLR 

IK 

DS  or  LS 

The  sehsnss  represented  in  Table  2  have  been  abbreviated  as  follows. 

IS  •  Zapliolt  lulsr  (Refs  17  and  14) 
n  •  Keller  box  (Ref  16) 

C8  •  Characteristic  box  <Rsf  19) 

08  ■  0wysr/8andsrs  <Refs  2  and  20) 

LS  ■  Laasonsn/Shsvslsv  (Refs  2,  21  and  22) 

LN  •  Lax/Wsndroff  type  (Refs  23  and  7) 

SB  •  Standard  box  (Ref  19) 

ZZ  •  Zig-zag  (Ref  24) 


4.2.3  Computational  grids  and  computing  effort 

In  a  coasHirison  of  boundary  layer  calculation  methods  carried  out  by  Eurovisc  at  the 
Amsterdam  Workshop*^,  the  required  computational  effort  was  deemed  an  important  indica¬ 
tor  of  program  performance.  As  computing  power  becomes  cheaper  this  is  a  measure  which 
is  receiving  less  emphssie  amongst  boundary  layer  method  developers,  but  nevertheless 
fcmsine  of  Interest,  not  leest  to  prospective  users.  The  speed  of  the  coaputer  in  appli¬ 
cations  depends  on  many  factors  amongst  which  may  be  mentioned  the  type  of  operation 
being  carried  out  and  how  the  program  being  run  is  organised.  The  charectsristlc  speed 
can  ba  axpraaaed  in  units  of  million  floating  point  operations  psr  second  (Mflope)  end 
is  evaluated  by  computing  standard  problems  with  the  same  teat  program  on  tha  different 
mmchinem  which  arc  to  bo  compared.  It  is,  however,  difficult  to  devise  a  test  program 
which  is  equally  fair  to  ell  of  the  various  machine  architectures.  In  othsr  words  thsre 
is  e  large  uneerteinty  (say  e  factor  of  two)  eeeocleted  with  the  resulting  cheracteri- 
etic  apaeda. 

Table  S  is  in  two  parte.  First  is  shown  the  coaputers  on  which  the  celculetlon 
methods  were  run,  together  with  estimated  charecterietic  machine  performance  figures  for 
eeientifie  eppiicatione.  The  second  pert  lists  the  execution  time  achieved  for  each 
repreeentetive  run  end  the  computational  grid  to  which  this  tins  relates.  Zn  the  lest 
column  Is  the  Effort  defined  ee 


Bffert 


Nf lope • tlme/nodes 


(4b:33) 
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Tabl*  3  ComPttftloiMil  effort 


Machine 

Speed 

(Hflops) 

Execution 
time  (see) 

Grid 

points 

Total 

nodes 

Effort 

X 

y 

z 

AVA 

IBM  4361 

0.86 

? 

100 

36 

o 

180*000 

7 

OAC 

CRAY  XMP/2 

24.0 

60 

7 

7 

o 

7 

7 

FFA 

Cyber  170/730 

0.22 

1790 

45 

37 

D 

80*000 

0.0048 

HBB 

VAX  8600 

0.84 

1900 

61 

81 

la 

335,000 

0.0038 

NLR 

Cyber  160/699 

2.0 

34 

40 

36 

O 

29,000 

0.0023 

Th«  flgur««  for  offort  aro  all  quit#  similar  but  tha  NLR  aathod  saama  to  stand  out  and 
It  Is  Indasd  trua  that  whan  It  was  daslgnad  a  dacada  ago  computational  afflclancy  was 
aada  an  Important  part  of  tha  opacification.  This  la  not  so  for  the  other  methods  in  tha 
table*  which  are  of  mors  recant  origin.  However*  It  Is  an  Interesting  Illustration  of 
Inherent  uncsrtalntlao  Involved  that  In  tha  Surovlsc  axarclaa^^.  referred  to  above*  the 
corresponding  Effort,  obtained  for  the  NLR  method  on  a  different  computer  was  found  to  be 
0.0092. 
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CHAPTER  5.  TEST  CASE  DATA  BASE 
5.1  Experimental 

5.1.1  Description  of  the  DFVLR  prolate  spheroid  and  its  instrumentation 

H.UJVleier  and  H.-P.Kreplin,  DLR,  Institut  fur  Exp.  Stromungsmechanik, 

BunsenstraBe  10,  D*3400  Gottingen,  Federal  Republic  of  Germany 


INTRODUCTION 

The  wind  tunnel  model,  especially  designed  for  the  investigation  of  three-dimensional 
boundary  layers  essentially  consists  of  a  1:6  prolate  spheroid.  This  configuration  was 
chosen  for  the  following  reasons: 

the  flow  pattern  around  an  inclined  prolate  spheroid  is  characteristic  for 
fuselages  and  missile  shapes, 

-  the  geometry  of  this  body  of  revolution  can  easily  be  described  applying 
orthogonal  or  non-orthogonal  coordinate  systems, 

-  the  potential  flow  is  given  -  even  for  the  inclined  version  -  by  closed  ana¬ 
lytical  solutions  [1], 

the  complete  instrumentation  like  probe  traversing  mechanism,  hot  film  ane¬ 
mometers,  pressure  transducers  etc.  can  be  storaged  in  the  model, 

as  the  model  can  be  turned  around  its  longitudinal  axis,  the  pressure  dis¬ 
tribution,  the  wall  shear  stress,  and  the  boundary  layer  quantities  can  be 
measured  in  each  cross  section  for  arbitrary  small  steps  of  the  circum¬ 
ferential  angle  $  with  the  same  surface  pressure  orifice,  surface  hot  film 
and  boundary  layer  probe,  respectively. 

1.  THE  WIND  TUNNEL  MODEL 

The  essential  construction  of  the  ellipsoid  is  based  on  two  half-shells  of  glass  fiber 
reinforced  resin  with  a  wall  thickness  of  about  4  mm.  The  dimensions  of  the  model  are 
given  in  Fig.  1,  which  shows  the  test  arrangement  in  the  3m  x  3m  Low  Speed  Wind  Tunnel 
of  the  DFVLR-AVA. 

The  model  can  be  turned  around  its  longitudinal  axis  through  <>  =  0*  to  360®  with  a  re¬ 
solution  of  A*  *  0.8®  by  means  of  a  step  motor.  For  one  particular  investigation  in  the 
DFVLR  3m  x  3m  Low  Speed  Wind  Tunnel  a  transition  strip  was  mounted  at  ?  0.2.  The 

strip  was  about  20  nun  wide  and  consisted  of  carborundum  (sand  paper  grain)  with  an  ave¬ 
rage  diameter  of  0.7  mm,  which  was  blown  onto  liquid  glue  on  the  model  surface  (for 
more  details  see  Ref.  C2] ) . 

The  construction  and  manufacturing  of  a  new  prolate  spheroid  for  the  investigations  in 
the  ONERA  FI  pressurised  wind  tunnel,  LeFauga,  became  necessary  due  to  the  seven  times 
higher  maximum  dynamic  pressure  in  the  FI  compared  to  the  Gbttingen  facility.  The  ob¬ 
jectives  of  the  DFVLR-ONERA  cooperation  were  to  extend  the  detailed  DFVLR  investigations 
on  the  prolate  spheroid  to  much  higher  Reynolds  numbers. 

2.  WALL  SHEAR  STRESS  MEASUREMENTS  BV  MEANS  OF  SURFACE  HOT  FILMS 

For  the  present  tests  two  prolate  spheroids  built  for  the  DFVLR  3m  x  3ir.  Low  Speed  Wind 
Tunnel  (NWG)  and  for  the  ONERA  FI  Wind  Tunnel  were  equipped  with  12  surface  hot  films. 

The  XQ-positions  of  the  flush  mounted  surface  hot  films  differ  only  slightly  on  the  two 
models  while  the  ♦-positions  are  completely  different,  see  Fig.  2.  The  commercially 
available  surface  hot  film  probes,  see  Fig.  3,  were  originally  designed  by  McCroskey  {3] 
for  the  Investigation  of  helicopter  rotor  blades.  Each  probe  consists  of  two  mutually 
perpendicular  films  in  V-conf iguration  on  a  thin  plastic  foil.  The  dimensions  of  the 
hot  films  are  given  in  Fig.  3.  Extreme  care  was  taken  to  mount  the  hot  films  on  the 
model  surface  in  order  to  avoid  an  imposed  boundary  layer  transition  caused  by  artificial 
roughness.  The  measured  imperfection  of  the  model  surface  in  the  regime  of  the  hot  films 
is  smaller  than  10  urn.  The  films  of  the  probes  were  operated  simultaneously  by  24  Con¬ 
stant  Temperature  Anemometers  developed  and  built  at  DFVLR,  Gbttingen.  These  anemometers 
were  located  inside  the  prolate  spheroid  in  order  to  avoid  very  long  probe  cables  in  the 
bridges.  The  overheat  ratio  of  the  films  was  adjusted  to  an  operating  film  temperature 
of  about  120®C. 

A  block  diagram  of  the  electronic  equipment  and  data  acquisition  system  used  for  the 
DFVLR  and  ONERA  tests  is  given  in  Figs.  4  and  5,  respectively.  A  detailed  description  of 
the  test  set  up  and  measuring  procedure  in  the  DFVLR  3m  x  3m  Low  Speed  Wind  Tunnel  is 
given  in  Refs.  [4],  C5]  and  for  the  ONERA  Fl  facility  in  Refs,  [6),  [71. 

Due  to  the  V-conf iguration  of  the  films  the  surface  hot  film  probes  can  be  used  for  the 
measurement  of  the  magnitude  and  direction  of  the  local  wall  shear  stress.  Changes  in 
the  local  wall  shear  stress  (magnitude  and/or  direction)  result  in  changes  of  the  local 
hot  film  transfer,  which  is  indicated  by  changes  of  the  anemometer  voltage  E.  The  heat 
transfer  per  unit  temperature  difference  is  used  in  the  calibration  and  data  reduction 
procedure.  This  quantity  is  given  by: 


Rp 

Q  =  = - =; -  with  E  =  I  R„  . 

‘f  ^Ref  * 

The  average  value  for  the  free  stream  and  wall  temperatures  was  used  as  the  reference 
temperature  T_  While  the  hot  film  resistance  R_  and  temperature  T  are  kept  constant 
by  the  anemometer,  the  reference  temperature  in  tne  wind  tunnel  may  change  with  time  so 
that  the  measured  free  stream  and  the  model  wall  temperature  correction  are  performed 
automatically,  which  is  sufficient  for  temperature  changes  of  a  few  degrees  as  observed 
in  the  experiments  reported  here. 

According  to  Ref,  C3]  the  wall  shear  stress  magnitude  is  derived  from  the  sum  of  the 
heat  transfer  rates  of  the  films  of  each  probe,  Q_  =  0^*^+  02*  The  measured  values  of  Qg 
were  related  to  calculated  wall  shear  stress  valuis  by'^the  following  procedure: 

For  zero  angle  of  incidence  the  mean  values  of  the  hot  film  signals  were  measured  at  dif¬ 
ferent  Reynolds  numbers  (Re  =  )•  These  measured  values  were  related  to  theoretical 

wall  shear  stress  values  i  ,  calculated  with  a  boundary  layer  computer  program  of  J.C. 
Rotta  [8]  for  the  given  fl<w  conditions.  The  calculations  were  based  on  the  measured 
pressure  distribution  and  the  experimentally  determined  transition  locations.  The  tran¬ 
sition  was  clearly  indicated  by  the  drastic  increase  of  the  anemometer  voltages  and  by 
fluctuations  in  the  hot  film  signals  without  a  calibration  of  the  probes.  Only  the  la¬ 
minar  and  fully  turbulent  boundary  layers  were  considered  for  the  calibration  procedure, 
because  the  transition  region  is  not  described  correctly  by  the  theory.  The  hot  film 
calibration  curves  obtained  in  this  axissymmetric  boundary  layer  flow  were  applied  for 
the  evaluation  of  the  local  wall  shear  stress  in  three-dimensional  boundary  layer  flows. 

In  order  to  be  able  to  take  small  variations  of  the  tunnel  temperature  into  account,  a 
calibration  of  the  surface  hot  films  was  performed  before  and  after  each  wall  shear 
stress  measurement  at  a  given  angle  of  incidence  and  free  stream  Reynolds  number. 

The  directional  sensitivity  of  the  hot  film  probes  could  not  be  derived  from  a  direct 
calibration  on  the  prolate  spheroid.  For  this  reason  the  calibration  was  carried  out  on 
a  flat  tunnel  wall  in  a  two-dimensional  turbulent  boundary  layer.  As  demonstrated  in 
Ref.  L 43  a  linear  relation  between  the  yaw  angle  t  and  the  hot  film  output  signals  was 
found  for  -40®  S  S  40®.  This  linear  relationship  was  applied  to  the  surface  hot  film 
probes  in  the  data'*  reduction  procedure. 

The  error  bound  for  the  magnitude  of  the  wall  shear  stress  depends  on  several  parameters 
like  Reynolds  number,  temperature  sensitivity  (substrate),  surface  curvature,  pressure 
gradient,  and  flow  direction.  For  the  reported  tests  the  estimated  accuracy  of  wall  shear 
stress  measurements  with  respect  to  the  magnitude  is 


Of  course,  the  variations  of  the  wall  shear  stess  can  be  measured  much  more  accurate 
(Within  the  order  of  1  %) .  The  determination  of  the  wall  shear  stress  direction  mainly 
depends  on  the  deviation  of  the  calibration  curve  from  its  linearity: 


ft  0.1  I Y  f  . 

w  ^  I  ’  w 


.  WALL  PRESSURE  MEASUREMENTS 

The  model  was  equipped  with  42  pressure  taps  positioned  on  one  meridian  in  non-equidis- 
tant  distances,  compare  Fig.  6.  Due  to  the  fact  that  the  model  could  be  rotated  around 
its  longitudinal  axis  in  steps  of  A*  i  0.8®  a  sufficient  spatial  resolution  in  circumfe¬ 
rential  direction  could  be  achieved. 

e  pressure  taps  have  a  small  diameter  (d  =  0.3  mm)  in  order  to  avoid  any  orifice 
effects  on  the  pressure  measurement  and  to  minimize  the  flow  disturbance  due  to  finite 
roughness.  The  pressure  taps  were  connected  to  a  pressure  scanner  (Scanivalve)  located 
inside  the  model.  For  the  experiments  carried  out  in  the  Gdttingen  wind  tunnel  (NWG) 
tne  integration  time  for  the  pressure  measurements  was  2  seconds.  The  pressure  inte- 
grat  ion  started  when  the  measured  pressure  value  was  indicated  to  be  constant  on  an  os- 
c.liosco[>e.  rt  was  checked  that  this  test  condition  was  reached  after  1  second.  This 
rtfeuJtfcd  ir,  a  necessary  time  for  the  measurement  for  all  wall  pressures  on  one  meridian 
•4rj  at  a  constant  circumferential  angle  of  about  2  minutes.  The  corresponding  pres¬ 
sures  were  measured  with  a  1000  Pa  pressure-transducer  manufactured  by  Hottinger 
(PUJ/0,M).  The  zero  shift  was  detected  before  each  run  =  const)  and  taken  into 
uccount  in  the  data  reduction  procedure.  Apart  from  small  temperature  changes  inside 
the  model,  the  zero  shift  of  the  pressure  transducer  was  mainly  caused  by  the  change 
of  ti.t  transducer  position,  if  the  model  was  rotated  or  inclined. 

hetails  about  the  test  arrangement  for  the  pressure  measurements  in  the  ONERA  FI  Wind 
iuj.nei  are  given  in  f6j. 

ihc  pressure  distributions  were  measured  in  the  Reynolds  number  range  of  Re  7  x  10^ 

(bf ’/I.R)  to  4Z  iO®  (ohrWA)  . 
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pressure  coefficient 


Cp  was  defined  as 


p  U^/2 


where  is  the  measured  wall  pressure  minus  the  static  pressure  and  p  U*/2  is  the 

dynamic  pressure,  respectively  denotes  free-stream  conditions).  The  static  pressure 
distribution  and  misalignment  of  the  free  jet  in  the  DPVLR  3m  x  3ra  Low  Speed  Wind  Tunnel 
NWG,  Gbttingen  is  given  in  Ref.  19).  The  measuring  accuracy  relative  to  Cp  =  1  is: 

flCp  =  0.005,  ONERA:  p^  =  3.85  bar,  =  75  m/s  and 

acp  =  0.01,  DFVLR:  =  1  bar,  *  55  m/s. 

4.  BOUNDARY  LAYER  MEASUREMENTS 

The  model  is  equipped  with  a  traversing  mechanism,  which  allows  vertical  and  horizontal 
movements  of  the  probe.  This  traversing  mechanism  is  installed  inside  the  model  in  such 
a  way  that  the  vertical  movement  relative  to  the  model  surface  is  warranted.  The  accu¬ 
racy  of  the  vertical  displacement  measurement  is  in  the  order  of  1/100  mm  and  the  hori¬ 
zontal  angle  measurement  is  better  than  0.4  degree. 

The  boundary  layer  measurements  were  carried  out  applying  a  Three-Hole-Direction-Probe 
with  an  outer  diameter  of  each  single  tube  of  d  =  0.5  mm  and  a  yaw  angle  of  45®,  com¬ 
pare  Ref.  [10).  This  probe  allows  the  determination  of  the  magnitude  and  the  direction 
of  the  local  velocity.  The  directional  sensitivity  with  respect  to  the  angle  of  inci¬ 
dence  a  and  the  yaw  angle  y  was  obtained  in  a  two-dimensional  channel  flow  of  the  DFVLR 
Low  Turbulence  Wind  Tunnel.  This  calibration  was  checked  during  the  tests  positioning 
the  probe  in  the  free  stream  while  the  angle  of  incidence  of  the  prolate  spheroid  was 
zero.  In  Ref.  [10)  the  directional  sensitivity  of  the  Three-Hole-Probe  for  yaw  angles  of 
y  -  ±  30®  and  angles  of  incidence  up  to  a  =  10®  is  given.  These  results  indicated  that 
the  directional  sensitivity  is  almost  independent  of  the  angle  of  incidence,  as  far  as 
they  are  representative  for  our  test  conditions.  This  is  not  true  for  the  total  pressure 
measurement  (p^) .  The  uncorrected  pressure  difference  p2“Pc»  non-dimensionalized  with 
the  dynamic  pressure,  obtained  at  free-stream  conditions,  changes  considerably  as  the 
probe  is  inclined  with  respect  to  the  tunnel  axis  to  a  *  10®.  However,  for  a  first  appro 
ximation  in  the  boundary  layer,  the  streamlines  are  parallel  to  the  model  surface  which 
implies  that  the  corrections  can  be  assumed  to  be  small.  Due  to  the  fact  that  during 
this  test  no  upwash  angles  were  measured  in  the  boundary  layers,  an  "a-correction"  was 
not  applied  in  the  data  reduction  procedure.  At  this  stage  of  investigation  the  static 
pressure  was  measured  at  the  model  surface  and  assumed  to  be  constant  through  the  en¬ 
tire  boundary  layer  thickness.  Errors  could  have  been  revealed,  here,  if  thick  boundary 
layers  close  to  separation  were  investigated. 

The  data  reduction  procedure  applied  is  described  in  Ref.  (10). 

Essentially,  the  resultant  velocity  can  be  derived  by  means  of  the  Bernoulli  equation, 
while  the  cross-flow  angle  y  is  given  by  the  pressure  difference  of  the  two  "45  degree 
tubes".  The  measurement  accuracy  of  the  resultant  velocity  relative  to  the  maximum  ve¬ 
locity  is 

fiu^  =  0.01 

while  the  directional  accuracy  is  better  than 

6y  -  I® 

in  the  linear  regime  of  the  calibration  curve.  For  wall  distances  of  z  <  1  mm  probe  dis¬ 
placement  effects  have  to  be  taken  into  account. 
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Fig.  1:  Prolate  spheroid  for  the  investigation  of  three-dimensional  boundary  layers 
in  the  3m  x  3m  Low  Speed  Wind  Tunnel  of  the  DFVLR,  G6ttingen. 


fig-  2; 


Length  of  the  model 

2a 

=  2.4 

m 

Diameter  of  the  model 

2b 

»  0,4 

m 

Maximum  angle  of  incidence 

a 

=  30® 

Maximum  velocity 

U^ 

=  62  m/s 

Turbulence  level 

TUj 

*  0.3 

% 

iQ[nvT\)  ift*) 
IH  0 

30  -IS 
S(7  .)0 

7S6  -tS 
962  0 

1166 
1370 
IS79 
1771 
1997 
2107 
2273 


-IS 

•30 


Locatic's  of  the  surface  hot  film  probes  in  different  cross 
prolate  spheroid  models  used  at  the  DFVLR  (a)  and  ONERA  (b) 
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Fi".  3:  Systematical  sketch  of  a  surface  hot  film  probe 


Computtr 

Ccfiltr 


Inltgrating 

Digital 

Vottmtttr 


24  Constant  • 
Tamptfotur*  > 
Anamomattfs 


24  Surfoca 
Hot  Films 


Fig.  4;  Block  diagram  of  the  electronic  equipment  and  data  acquisition  system 
used  in  the  DFVLR  3m  x  3m  Low  Speed  Wind  Tunnel 
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Fig,  5;  Block  diagram  of  the  electronic  equipment  and  data  acquisition  system 
used  in  the  ONERA  FI  Wind  Tunnel 
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5,1.2  Lockheed/NASA-Ames  Wing  C 

H,  Yoshihara 
Boeing  Company 
Seattle,  Washington 


The  second  case  chosen  by  VG  10  is  tfing  C  at  H  •  0,85,  5  degrees  angle  of  attack,  and  Re,  no,  -  6.8  x 
10^/aean  chord.  Vind  tunnel  aeasureaents  of  the  surface  pressures  and  boundary  layer  velocity  profiles 
have  been  obtained  by  Keener  (Refs.  1  and  2).  In  the  folloving  the  definition  of  the  ving,  details  of  the 
test,  and  the  plots  and  tabulation  of  the  test  data  taken  froa  Refs.  1  and  2  are  presented.  In  addition, 
calculations  vith  the  Transonic  Navier/Stokes  (TNS)  code  (Ref.  3)  are  given  to  illustrate  the  possible 
test/theory  eisaatches  that  can  be  expected. 

2.1  Details  of  the  Test  and  Test  Data 

The  tests  vere  carried  out  in  the  Aaes  6x6  Foot  tunnel  vhich  is  a  continuous  flov  facility  having 
slotted  upper  and  lover  surface  vails  of  6X  porosity,  Geoaetric  properties  of  the  aodel  are  suaaarized  in 
Figure  1.  The  ving  vas  defined  in  teras  of  the  root  and  tip  chords  vith  the  interaediate  sections  defined 
by  a  linear  developaent.  The  root  and  tip  sections  are  shovn  in  Figure  2,  and  their  coordinates  are  given 
in  Table  1.  Since  a  definition  of  the  ving  tip  vas  not  provided,  a  shape  as  shovn  in  Figure  1  vas 
postulated.  The  ving  vas  aounted  on  the  tunnel  aidevall  directly  as  a  seai^span  aodel.  Its  seal-span  of 
0.9  aeters  resulted  in  a  test-section  blockage  ratio  of  1.3X.  Details  of  the  construction  tolerances  and 
the  aeroelastic  properties  are  given  in  Ref.  1.  Surface  pressures  vere  obtained  by  a  transducer- 
scanivalve  systea,  ^ile  the  boundary  layer  velocity  profiles  vere  obtained  by  a  3-hole  "cobra  head”  flov 
direction  probe  vith  the  traversing  aechanisa  aounted  through  the  ving.  The  location  of  the  surface 
pressure  taps  and  the  velocity  profile  stations  are  shovn  in  Figure  3.  The  boundary  layer  vas  tripped 
using  sifted  glass  spherules  located  at  the  4,S!t  chord  station.  The  final  deteraination  of  the  trip 
configuration  vas  deterained  using  subliaation  tests.  On  the  lover  surface  and  outboard  of  the  602  span 
station  on  the  upper  surface  the  trip  size  vas  0.16  an,  vhile  inboard  of  the  602  span  station  on  the  upper 
surface  the  trip  size  vas  0.23  aa.  Fluorescent  oil  flov  pictures  vere  also  obtained. 

The  test  data  for  K  >  0.85  and  5^  angle  of  attack  are  plotted  in  Figures  4  and  5  and  tabulated  in  Tables  2 
and  3. 

2.2  Transonic  Navier/Stokes  Solution 

Calculation  for  the  above  case  using  the  TNS  code  vas  carried  out  earlier  (Ref.  3).  The  Baldvin/Loaax 
equllibriuB  turbulence  aodel  vas  eaployed,  and  an  approxiaate  syeaetry  condition  vas  iaposed  at  the  root 
chord  plane.  The  calculations  vere  therefore  repeated  vith  the  exact  syaaetry  condition  iapleaented  using 
an  iaage  plane.  In  Figure  6  the  resulting  chordvise  pressure  distributions  are  coapared  to  the  earlier 
Ref.  3  result  and  vith  the  experiaental  result.  Coaparison  of  the  tvo  calculated  results  shovs  that  the 
suctions  are  significantly  higher  vith  the  exact  symaetry  conditions  in  the  upper  surface  peak  and  plateau 
regions  vith  a  smaller  increase  of  the  suction  plateau  on  the  lover  surface.  These  effects  of  the 
syaaetry  condition  extend  over  the  entire  span,  such  tac-ceaching  Influences  being  characteiistic  of  high 
subsonic  flovs. 

Surprisingly  the  less  exact  calculations  agreed  aore  closely  vith  the  experlaents  so  far  as  the  suction 
peaks  and  plateaus  vere  concerned,  A  partial  explanation  for  this  is  the  distortion  of  the  test  data  by 
the  significant  boundary  layer  displaceaent  effect  on  the  sidevall  on  vhich  the  ving  is  aounted.  Such 
sidevall  displaceaents  vill  be  coapressive  about  the  ving  leading  edge,  depressing  the  suctions  in  the 
upper  surface  peak  and  plateau  regions  over  the  entire  span  in  auch  the  saae  Banner  as  the  inexact 
syaaetry  conditions.  In  addition  the  consequences  of  the  siaplifications  used  in  the  diagonalization 
process  in  the  TNS  code,  particularly  in  the  leading  edge  region,  need  to  be  assessed.  The  remaining 
large  test/theory  alsaatches,  as  at  the  702  and  902  span  stations,  are  most  probably  due  to  the  inadequacy 
of  the  turbulence  aodel  and  to  the  inadequate  treatment  of  the  tip  flov.  Thus  for  example  the  difference 
of  the  pressure  distribution  at  the  702  station,  namely  the  calculated  single  shock  distribution  compared 
to  the  measured  double-shock  distribution,  stems  froa  the  inboard  displacement  of  the  forvard  and  rear 
shock  intersection  by  the  inadequate  turbulence  aodel.  The  aatter  of  the  turbulence  modeling,  however, 
falls  outside  the  scope  of  the  VG  10. 

The  above  calculations  vere  carried  out  by  G.  Blom  and  T.  Blum  of  the  Boeing  Coapany. 
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1.  Introduction 

The  prolate  spheroid  is  a  body  of  simple  geometry  that  can  provide  significant  understanding  of  the  complex  flow 
that  can  develop  about  slender  bodies  at  incidence.  A  variety  of  wind  tunnel  tests  have  been  carried  out  for  a  0:1  prolate 
spheroid  at  the  DFVLR,  at  Gottingen,  using  a  glass-fiber  model  of  4mm  wall  thickness  and  having  major  and  minor  axis 
of  2.4m  and  0.4m,  respectively.  The  majority  of  the  tests  have  been  carried  out  at  moderate  Rejmolds  numbers  but  some 
measurements  have  also  been  conducted  in  collaboration  with  ONERA,  in  FVance,  at  higher  Reynolds  numbers. 

For  the  purposes  of  the  Working  Group  No  10  of  the  Fluid  Dynamics  Panel  of  AGARD,  the  flow  field  surrounding 
a  fi:l  prolate  spheroid  was  calculated  at  the  NASA  Ames  Research  Center.  Two  sets  of  flow  conditions  were  studied, 
a  =  30**,  (/oo  =  75m/6,  Rej  =  43x10*,  corresponding  to  conditions  tested  at  ONERA,  and  the  transitional  flow  case  of 
a  =  10”,  Ugo  =  SSm/s  and  Aei  =  7.7x10*,  which  the  members  of  the  Working  Group  have  used  for  the  boundary-layer 
calculations.  For  the  calculation  of  the  flow  field  the  code  F3D  described  in  Ref.  1  has  been  used.  The  code  solves  a 
conservative  form  of  the  thin-layer  Navier-Stokes  equations  cast  in  ^neralised  coordinates,  so  that  it  can  be  readily  used 
for  computing  flows  about  complex  configurations. 

Here  the  results  of  these  calculations  and  comparisons  with  the  available  experimental  data  are  pven.  The  calculated 
flow  fields  are  similar  to  the  ones  obtained  experimentally,  characterised  by  the  appearance  of  a  pair  of  symmetrical  vortices 
at  the  leeward  of  the  body,  that  ori^nates  along  the  lines  of  cross-flow  separation.  Similar  computational  results  for  these 
cases  were  also  recently  published  by  Vatsa  et  al.  (Ref.  2). 

2.  Numerical  method 

The  conservation  equations  of  mass,  momentum,  and  energy  can  be  represented  in  a  flux- vector  form  that  is  convenient 
for  numerical  simulation  as: 

drQ  +  &((F  +  fV)  +  d,(5  +  Gg)  +  d((ff  +  5,)  =  0  (1) 

where  r  is  the  time,  and  the  independent  spatial  variables  ^,17,  and  (  are  chosen  to  map  a  curvilinear  body- conforming 
discretization  into  a  uniform  computational  space.  In  the  system  used  in  this  study,  (  denotes  the  curvilinear  axis  in  the 
direction  of  the  main  l^dy  sods,  rf  denotes  the  circumferential  and  (  the  normal  to  the  body  direction.  As  opposed  to  the 
invisdd  flux  terms  F,G,  and  B,  the  terms  Fp^Gg  and  Hp  are  fluxes  containing  the  viscous  derivatives.  A  nondimensional 
form  of  the  equations  is  used  throughout  this  work.  Lengths  are  scaled  by  the  length  of  the  body,  L,  velocity  components 
by  the  free-stream  velocity  of  sound,  Ooq,  the  pressure  and  the  total  energy  per  unit  volume,  (e),  by  paeOoo^*  other 
quantities,  such  as  T,  p,  m  are  scaled  by  their  &ee-stream  values. 


For  body-conforming  coordinates  and  high-Reynolds  number  flow  the  thin-layer  approximation  can  be  applied  (Ref. 
3,4) 

arQ  +  a(F  +  e,d  +  d(B  =  iu-'d(S  (2) 

where  only  viscous  terms  in  the  normal  to  the  body  direction  are  retained.  These  have  been  collected  into  the  vector  S  and 
the  nondimensional  Reynolds  number  Re  is  factored  from  the  viscous  flux  term.  In  the  present  turbulent  computations, 
the  coeflSdents  of  viscosity  appearing  in  the  right-hand-side  terms  of  Eq.  (2)  are  obtained  from  the  modd  of  Baldwin  and 
Lomax  (Ref.  3). 

In  differendng  these  equations  it  is  advantageous  to  difference  about  a  known  base  solution.  If  the  base  state  is 
properly  chosen,  the  differenced  quantities  can  have  smaller  and  smoother  variation  and  therefore  less  differendng  error. 
If  the  free-stream  is  used  as  the  base  solution: 

+  (((H-S„)-Re-'S((§-S^)  =  0 

where  6  indicates  a  general  difference  operator. 

An  implidt  approximately-factored  scheme  for  the  thin-layer,  Navier-Stokes  equations  that  uses  central  differendng 
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ia  the  17  &nd  ^  directions  and  upwind  differencing  in  the  ^  direction  can  be  written  in  the  form: 

[/  +  +  htfC’  -  -  Di\(] 

X  (/  +  hS{(A- )"  +  -  DilJAQ*  = 

-  A<{f}[(F+)"  -  /^] + </i(F-r  -  ^i;;i  w 

+  i,(G"  -  Goo)  +  S^S'  -  5oo)  -  lU-'idS^  -  Soo)} 

-D.($"-Qoo) 

where  h  =  At  or  At/2  for  first  or  second  order  time  accuracy,  and  the  free  stream  base  solution  is  used.  Here  6  is  typically 
a  three-point,  second-order-accurate,  central  difference  operator,  5  is  a  midpoint  operator  used  with  the  viscous  terms, 
and  the  operators  and  6^  are  backward  and  forward  three-point  difference  operators.  The  fiuz  F  has  been  split  into 
and  F~,  according  to  its  eigenvalues  (Kef.  5),  and  the  matrices  and  Af  result  &om  local  linearisation  of  the 

fluxes  about  the  previous  time  level.  In  Eq.  (4),  J  denotes  the  Jacobian  of  the  coordinate  transformation.  Dissipation 
operators,  D*  and  Dj,  are  used  in  the  central  space  differencing  directions  (Ref.  1). 

The  factored  left-hand-side  operators  can  be  readily  solved  by  sweeping  in  the  i  direction  and  inverting  tridiagonal 
matrices  with  5x5  blocks.  This  two-factor  implicit  scheme  is  readily  vectorised  in  planes  of  ^  r=:  constant. 

The  grid  was  generated  by  solving  a  system  of  hyperbolic  differential  equations  (Sieger  and  Chaussee,  Ref.  6).  First, 
a  two-dimensional  grid  was  generated  in  the  symmetry  plane.  Then  the  three-dimensional  grid  was  obtained  by  rotating 
it  about  the  major  axis  of  the  ellipse.  The  grid  consists  of  121  points  along  the  direction  of  the  main  axis  of  the  body, 
100  around  its  circular  section  (no  plane  of  63rmmetry)  and  65  radially  from  the  surface.  For  the  various  graphics,  the 
output  file  was  reduced.  The  reduced  grid  is  shown  in  figure  1.  grid  point  (J,  K,  L)  is  defined  by  the  Cartesian 

coordinates  x,  y,  s.  The  axial  coordinate,  x,  is  measured  from  the  nose  of  the  ellipsoid,  the  lateral  coordinate,  y,  from  the 
vertical-longitudinal  symmetry  plane  and  the  normal  coordinate,  s,  from  the  horizontal-longitudinal  symmetry  plane. 


3.  Calculated  quantities 

For  the  calculation  of  the  skin  friction  the  physical  contravariant  velocity  components  ,U*,  V*,  W*,  were  used. 
The  physical  contravariant  velocity  components  are  related  to  the  cartesian  velocity  components  through  the  following 
equations: 

V"  =  (u{,  +  vi,  + 

V’  =  (ui;,  +  tiT),  +  +  + 

W'  =  (u<.  +  vi,  +  (5) 

The  vectors  etc.,  are  the  components  of  the  contravariant  base  vectors,  while  the  vectors  etc.,  are 

the  components  of  the  covariant  metrics.  Equations  (5)  are  valid  for  a  curvilinear  non-orthogonal  system. 


The  components  of  the  local  skin-friction  coefficient,  along  the  streamwise  and  along  the  crosswise  directions,  were 
estimated  by  the  equations: 

"  900^'' 

~  ac 


As  the  skin  friction  is  estimated  near  the  wall,  only  the  djmamic  viscosity  fi  was  considered.  Its  value  was  estimated 
by  Sutherland's  formula.  The  velocity  differences  and  the  associated  normal  distance  were  estimated  between  the  wall  (L 
s  1)  and  the  first  layer  of  grid  points  (L  —  2). 


4.  Results 

Morphology  of  the  flow  surrounding  inclined  bodies 

It  is  known  that  the  flow  field  surrounding  bodies  at  incidence  is  characterised  by  the  appearance  of  a  leeward-side 
vortex  structure.  The  morphology  of  the  vortical  structure  is  similar  for  a  variaty  of  shapes  of  bodies.  This  topic  has  been 
covered  in  detail  by  Tobak  and  Peake  (Ref.  7)  and  by  Perry,  Chong  and  Bomung  (Ref.  8).  In  a  cross-section  of  such  a 
flow,  under  the  action  of  circumferential  pressure  gradient,  the  outer  flow  approaching  the  windward  plane  of  sjrmmetry 
toms  and  flows  outwards  along  the  body,  from  the  windward  toward  the  leeward  side.  The  boundary  layer  which  it 
formed  in  this  way,  separates  from  the  body  at  a  point  on  the  leeward  aide  that  depends  on  the  streamwise  position  of 
the  cross-sectiOT  and  on  the  flow  condiiioni.  Then  the  fluid  leaves  the  body  along  a  feeding  sheet  and  rolls  up  to  form  a 
primary  vortex  system  on  the  leeward  side  of  the  body.  The  pair  of  the  primary  vortices  induces  a  flow  toward  the  body 
surface,  which  at  the  point  of  attachment  turns  outwards  toward  the  windward  side.  For  primary  vortices  of  sufficient 
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•treogthf  a  secondary  separation  is  indnced  to  the  boundary  layer  below  them.  Thus  the  secondary  Tortex  structure  may 
appear,  that  rotates  in  an  opposite  direction  than  the  primary  one.  Both  the  systems  are  carried  donwstream  by  the  axial 
component  of  the  flow. 

Experimentally,  the  vortex  structures  may  be  detected  by  particle  trace  techniques  (e.g.  smoke),  reconstruction  of 
velocity  measurements,  or  more  simply,  by  using  sur&ce  <^>flow  techniques.  The  lines  that  are  formed  on  a  surface  by  the 
oil  correspond  to  the  surface  sJdU'&iction  lines.  These  lines  converge  toward  the  primary  and  the  secondary  separation  lines 
and  diverge  from  the  attachment  line,  that  lies  between  them.  Computationally,  the  equivalent  of  the  surface  shear-stress 
lines  visualised  in  oil-flow  experiments  is  simulated  by  releasing  particles  at  selected  grid  points  on  the  first  layer  of  pmnts 
above  the  body.  The  particles  move  with  the  local  fluid  velocity,  but  are  constrained  to  the  plane  just  above  the  body,  i.e, 
the  normal  velocity  component  is  set  to  sero  in  the  time  integration. 

Comparison  of  the  results  with  experiments,  a  =  30^ 

The  measurements  at  ONERA  included  surface- pressure  and  skin  friction  distribution,  together  with  oil-flow  visualisa¬ 
tion.  The  more  detailed  tests  performed  at  Gottingen  at  a  Reynolds  number  equal  to  8.52x10*  also  include  measurements 
of  the  velocity  components  in  the  flow  field  (Ref.  9).  The  experimental  results  indicated  that  in  both  cases  the  flow 
developed  similarly  to  the  typical  case  described  in  the  previous  section.  At  this  time  the  ONERA  experimental  data  have 
not  been  published.  However,  Meier  (private  communication)  has  provided  the  distribution  of  the  pressure  coefficient  and 
some  other  information. 

Flow-field  results  of  the  calculations  are  shown  in  figures  2,  3.  The  existence  of  the  vortex  structure  and  its  outer 
extent  at  midlength  of  the  body  is  indicated  by  the  velocity-vectors  plot  (figure  2).  The  primary  and  secondary  vortex 
structures  are  detected  in  figure  3,  where  surface  particle  traces  are  used  to  rimulate  the  surface  shear-stress  lines  visualixed 
in  oU-fiow  experiments.  In  this  figure  the  experimental  oil-flow  pictures  included  in  Ref.  10  are  also  shown,  in  top  and  side 
view.  It  is  observed  that  the  computed  particle  traces  show  relatively  good  agreement  with  the  experiments.  The  extent 
of  the  crossflow  separation,  indicated  by  the  primary  separation  line,  is  very  similar  in  both  cases,  with  the  exception  of  a 
small  re^on  near  the  sting.  There  the  separation  is  more  extensive  for  the  calculated  results  (inadvertent  use  of  a  thinner 
sting  in  the  computations  may  account  for  some  of  this  discrepancy).  Also,  the  coincidence  of  the  secondary  crossflow 
separation  lines  is  remarkable.  In  both  casCb  they  run  along  the  azimuthal  line  of  30**  from  the  leeward  symmetry  plane. 

Computed  surface-pressure  coefficients  at  windward  (^  =  0^ )  and  leeward  (^  s  ISO**)  planes  of  symmetry  are  compared 
with  the  experimental  data  in  figure  4.  The  agreement  with  the  experimental  results  on  the  windward  side,  it  good.  This 
was  expected,  because  the  windward  boundary  layer  remains  attached  almost  over  the  entire  length  (near  the  sting  there 
are  again  some  differences).  The  computed  and  experimental  pressure  distributions  are  also  in  good  agreement  on  the 
leeward  plane  of  symmetry.  A  small  difference  is  observed  between  x/L  =  0.2  -  0.4.  It  is  noted  that  while  in  the  experiments 
the  flow  was  transitional  at  part  of  the  nose,  io  the  calculations  it  has  b<^  assumed  that  the  flow  starts  as  a  turbulent 
one. 


In  order  to  verify  the  quality  of  the  numerical  solution  away  from  the  symmetry  planes,  the  computed  and  the 
experimental  pressure  coefficients  are  compared  along  the  drcumferential  direction,  in  figure  5,  at  three  axial  locations. 
The  agreement  is  good,  not  only  on  the  windward  but  also  on  the  leeward  side.  The  suction  observed  at  ^  160*,  on  the 

two  upper  curves,  is  due  to  the  secondary  vortex. 

The  computed  shear  stress  angle,  7,  along  the  wall,  is  shown  in  figure  6,  at  the  same  axial  locations  that  were  used  in 
the  case  of  the  pressure  coeffident.  The  experimental  data  have  not  been  released  as  yet,  so  a  comparison  is  not  possible. 
However  the  good  agreement  of  the  oil-flow  simulation  with  the  experiments  indicates  that  the  predicted  7  angles  may  be 
not  very  different  from  the  experimentally  found  values. 

Effect  of  turbulence  model  and  of  grid  refinement 

Both  the  turbulence  modeling  and  grid  resolution  can  affect  the  output  of  the  calculations  and  must  be  assessed.  In 
the  present  paper,  the  calculations  were  started  with  a  gnd  consisting  of  93x72x47  points  and  then  they  continued  with  a 
grid  of  121x100x65  points. 

I”  the  case  of  the  turbulent  modeling,  the  tiirbulent  viscosity  coeffident  ih  was  computed  using  the  two-layer,  Cebed- 
type,  .ugebraic  eddy-viscosity  model  reported  by  Baldwin  and  Lomax  (Ref.  3).  In  both  layers,  depends  on  the  absolute 
value  of  the  local  vortidty  vector,  uf,  and  on  the  distance  from  the  wall,  y.  More  spedfically,  in  the  outer  layer  the  turbulent 
viscosity  coeffident  is  proportional  to  the  maximum  value  of  the  function: 

F(v)  =  |«|vll-e-<J^>|  (7) 

and  io  the  value  of  y  at  which  this  maximum  occurs,  ym««< 

Degani  and  Schiff  (Ref.  11)  observed  that  the  evaluation  of  yms*  i«  not  straightforward  in  cases  of  bodies  with 
crossflow  separation,  bideed,  while  in  the  windward  dde  the  attached  boundary  layer  gives  rise  to  a  pr^e  of  F{y)  which 
has  a  tingle  maidmum,  in  the  leeward  side  a  second  nMudmnm,  of  greater  value,  appears.  This  second  maiitnnm  is  due  to 
the  overiying  vortical  structure.  Thus,  if  in  the  code,  the  computer  searches  outward  along  each  ray  to  determine  the  peak 
of  F'(y),  it  will  select  the  second  maximam.  The  resulting  value  of  will  be  much  higher  than  the  one  of  the  boundary 
layer.  This  will  cause  an  underestimation  of  the  extent  of  the  crossflow  separation  and,  consequently,  of  the  die  of  the 
vortical  structures. 
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For  eliminating  the  forementioned  difficulty,  Degani  and  Schiff  modified  the  turbulence  subroutine,  in  their  code,  so 
that  the  first  peak  of  F(y)  is  selected  at  each  profile.  In  the  present  calculations,  initially  the  flow  was  estinoated  without 
the  modification  suggested  by  Degani  and  Schiff.  Then  after  reviewing  the  profiles  of  the  function  F'(y)  around  the  prolate 
spheroid,  a  cut-off  distance  was  defined,  that  divides  the  bounda^  layer  and  the  vortical  structures.  This  is  evident  in 
figure  7,  where  profiles  of  F(y)  are  shown  at  some  typical  points  of  the  body.  It  is  observed  in  this  figure  that  the  maxima 
of  the  boundary-layer  vortidty  distributions  lie  below  the  value  L  =  22,  in  all  the  stations.  Thus,  if  the  outward  search 
stops  at  this  cut-off  distance  only  the  first  peak  will  be  selected  and  not  the  one  which  is  due  to  the  vortical  structures. 
After  this  principle  was  applied  for  the  evaluation  of  the  turbulent  viscosity  coefficient,  the  profiles  of  F(y)  were  plotted 
again,  and  checked  as  to  whether  the  value  L  =  22  was  stiU  appropriate.  No  significant  change  has  been  noticed. 

The  above  observation  helps,  also,  in  defining  the  edge  of  the  boundary  layer.  It  is  evident,  in  figure  7,  that  the 
voxiidty  function  F(v),  after  it  takes  its  local  maximum  value  at  the  boundary  layer,  it  falls  rapidly  towards  a  minimum 
value  and  then  it  starts  again  to  rise,  at  the  leeward  side,  where  the  vorti^  are  present.  At  the  windward  side,  where  the 
boundary  layer  is  attached,  only  one  maximum  appears  and  then  the  F(y)  falls  to  zero.  So  in  general  the  function  F(y) 
indicates  rather  dearly  the  position  of  the  edge  of  the  boundary  layer. 

Although  a  systematic  study  of  grid-refinement  and  variation  of  the  turbulence  model  was  not  carried  out,  computed 
surface  streamline  patterns  are  shown  in  figure  8,  for  a  coarse  grid(8a),  a  fine  one  with  the  standard  turbulence  model(8b), 
and  a  fine  one  with  the  improved  turbulence  model(8c).  For  a  better  comparison,  an  tmwrapped  coordinate  system  is 
used.  Note  that  in  this  land  of  graphics  the  pictures  are  distorted,  because  the  vertical  absissa  takes  the  same  value  at 
the  position  of  the  maximum  diameter  and  at  the  position  of  the  sting  or  the  nose.  A  comparison  shows  that  as  far  as 
the  primary  separation  line  is  concerned,  the  three  calculations  show  essentially  the  same  results.  However,  things  are 
completely  different  in  the  case  of  the  secondary  separation  line.  In  the  coarse-grid  calculation,  this  line  does  not  appear 
at  ail.  In  the  fine-grid  calculations  with  the  unmodified  turbulence  model,  an  attachment  region  has  been  formed  between 
the  primary  separation  line  and  the  symmetry  plane.  In  the  final  calculations,  after  the  improvement  of  the  turbulence 
model,  the  secondary  separation  line  appears  very  clearly. 

Finally,  the  circumferential  variation  of  the  pressure  coefficient  at  the  location  x/L  =  0.81  is  shown  in  figure  9  for  the 
three  calculations.  The  experimental  data  are  also  included.  The  gradual  improvement  of  the  predicted  results  is  evident 
in  this  figure,  also. 


Compairison  of  the  results  with  experiments,  a  =  10* 

The  experiments  for  the  a  =  10*  case  were  conducted  at  the  DFVLR.  The  Reynolds  number  based  on  the  major  axis 
of  the  ellipsoid  was  7.7  x  10*.  The  tunnel  velocity  was  fiSm/s.  Transition  was  imposed  at  20%  of  the  major  axis  from 
the  nose,  by  a  circumferential  strip  of  carborundum  powder.  The  tests  indicated  a  fully  turbulent  boundary  layer  next  to 
the  transition  strip.  Meier  (ref.  10)  has  given  the  experimentally  obtained  surface  shear-stress  and  flow  inclination  data. 
Velocity  profiles  at  various  stations  are  included  in  Ref.  12. 

For  the  calculation  it  was  assumed  that  the  flow  was  laminar  up  to  x/L  =  0.20  At  this  point  an  instantaneous  transition 
to  turbulent  flow  was  imposed.  Surface  particle-traces  of  the  computed  solution  are  shown  in  figure  10.  It  is  observed 
that  the  cross-flow  separation  is  limited,  compared  to  the  30*-case.  Also  no  secondary  vortex  structure  appears,  evidently 
because  the  strength  of  the  primary  vortex  is  not  sufficiently  strong  to  induce  it.  The  experimental  oil-flow  results  are  not 
available.  However,  the  crossflow  angle,  7,  is  known.  The  calculated  crossflow  angle  and  the  experimentally  found  values 
are  shown  in  figure  11,  at  various  cross-sections.  Like  the  pressure  distribution,  the  crossflow  angle,  7,  is  much  easier  to 
accurately  predict  than  the  skin  friction. 

The  computed  sldn-friction  coefficient  distribution,  in  the  circumferential  direction  at  three  axial  stations,  is  compared 
with  the  experimental  data  in  figure  12a.  The  range  of  the  parameter  of  the  flow,  on  the  first  layer  of  points  above  the 
surface  of  the  body  was  y'*'  =  0.5-0.9.  The  agreement  of  the  predicted  data  with  the  experimental  values  is  quite  good. 
Shown  in  figure  12b  is  a  general  coordinate  boundary  layer  result  (Ref.  13)  ,  in  which  edge  conditions  were  taken  from  the 
Navier  Stokes  result  at  approximately  30mm  above  the  wall.  The  boundary  layer  code  uses  central  differencing  without 
added  numerical  dissipation  in  the  direction  away  from  the  wall,  and  it  is  used  to  verify  the  Navier-Stokes  predicted  values 
of  skin  friction.  Previously  obtained  Navier-Stokes  results  using  more  numerical  smoothing  and  a  coarser  grid  produced 
different  level  of  C/,  while  the  boundary  layer  code  was  not  sensitive  to  those  changes,  even  though  it  used  the  different 
Navier-Stokes  inputs. 

Velocity  profiles,  on  the  symmetry  plane  at  x/L  =  0.65  and  0.74,  are  compared  with  the  experiments  in  figure  13a. 
The  agreement  on  the  windward  side  is  better  than  00  the  leeward  side,  where  the  measured  profiles  are  fuller  in  their 
lower  part  than  the  predicted  ones.  For  comparison,  the  computed  boundary  layer  profiles  01  Ref.  13  are  shown  in  figure 
13b. 


In  first-order-of-magnitude  boundary-layer  methods  the  assumption  it  made  that  the  pressure  is  constant  across  a 
boundary  layer.  To  check  this  assumption  for  the  particular  prolate  spheroid  flows  that  are  reported  in  the  present  pnper, 
the  wall  pressure  in  the  drcumferential  direction  (x/L  —  0.74)  is  compared  in  figure  14  with  the  preMure  at  a  normal 
distance  from  the  wall.  This  distance  is  of  the  order  of  the  maximum  thickness  of  the  boundary  layer  (30mm).  It  is  teen, 
in  figure  14,  that  the  pressure  is  nearly  constant  irithin  the  boundary  layer,  in  both  the  10*  and  30*  case.  A  small  gradient 
^>pears  only  in  the  30*  case  at  the  region  of  maximum  accelerati<m  of  the  croN  flow. 
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Fig.  1.  Schematic  view  of  the  grid.  Only  a  few  lines  are  shown. 
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Fig.  2.  The  calculated  flow  field.  Crossflow  velocity  vectors  at  x/L 
Rei  =  43x10®,  a  =  30°. 


Fig.  5.  Compaiison  of  surface-pressure  at  various  cross-flow  sections. 
Re,  =  43x10*,  a  =  30". 
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Fig.  7.  Profiles  of  the  vortidty  function,  F(y).  (a)  axial  location 
x/L=  0.35.  (b)  axial  location  x/L=  0.90. 

Rei  =  43xl0»,  a  =  30'. 
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Fig.  9.  Comparison  of  the  calculated  surface  pressure  in  the  circum¬ 
ferential  direction  at  x/L=  0.81. 

Rei  =  43x10“,  a  =  30". 
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Fig.  10.  Surface  paiticle-tiaces.  Rei  =  7.7x10*,  a  =  10'. 
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Fig.  13.  Velocity  profiles  comparison  at  the  symmetry  planes,  (a) 
x/L=0.65.  (b)x/L=0.74. 

Rei  =  7.7xl0»,  a  =  10”. 
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CHAPTER  6 .  RESULTS 

To  investigate  the  accucacy  o£  pcedicting  thcee-dimensional  flows  by  solving  the 
boundary-layer  equations  with  field  and  integral  sethods,  two  test  cases  have  been  con¬ 
sidered  as  a  result  of  the  discussions  of  Chapter  S.  The  first  corresponds  to  the  pro¬ 
late  spheroid  discussed  in  Chapter  6  and  exaained  experimentally  In  the  3  m  x  3  m  wind 
tunnel  of  the  DPVLR  in  Gdttlngen.  Measurements  of  wall  shear  and  velocity  profiles  were 
obtained  at  four  streamwise  stations.  x/2a  «  0.4S.  0.S6.  0.64  and  0.73  for  an  angle  of 
attack  of  o  >  10*.  The  measured  distributions  of  the  wall  pressure  coefficient  provide 
a  basis  for  determining  the  external  velocity  boundary  condition  required  in  the  solution 
of  the  boundary- layer  equations.  An  alternative  approach  to  this  boundary  condition  is 
to  make  use  of  the  measured  external  velocity  values  but  these  are  not  provided  in  suf¬ 
ficient  detail  so  that  other  means  require  investigation. 

The  uncertainties  associated  with  the  specification  of  the  external  velocity  dis¬ 
tribution  have  been  investigated  with  results  shown  in  Pig.  I.  These  results  correspond 
to  the  velocity  components  u^/u*  and  We/u«  at  x/2a-8tatione  of  0.48  and  0.73  obtained  by 
the  solution  of  the  inviscid  flow  equations  with  the  surface  pressure  measurements  com¬ 
puted  by  Stock,  and  from  the  measurements  of  Meier  and  Kreplin.  It  is  clear  that  the 
results  of  the  inviscid  flow  theory  are  usually  higher  than  those  computed  from  the  wall 
pressure  distribution  and  considerably  higher  than  those  measured  at  the  edge  of  the 
boundary  layer.  One  possibility  for  these  discrepancies  may  be  the  pressure  variation 
across  the  boundary  layer,  but  this  cannot  be  entirely  true  since  the  difference  in  vel¬ 
ocities  also  appears  in  areas  where  the  boundary  layer  is  thin  and  laminar,  for  example, 
near  the  nose.  Another  possibility  may  be  that  the  fiowfield  in  the  open  test  section 
is  deflected  by  the  presence  of  the  body  because  it  imposes  lift  and  consequently  dis¬ 
torts  the  pressure  distribution.  In  the  light  of  these  uncertainties,  the  Working  Group 
defined  a  standard  test  case,  at  10*  incidence,  which  used  the  external  flow  data  as 
derived  from  the  wall  pressure  distribution  by  Stock. 

In  general  the  velocities  calculated  from  the  wall  pressure  distribution  by  Stock 
are  lower  than  those  calculated  from  the  potential  flow  but  the  difference  is  only  a  few 
percent  for  4  <  120*.  For  4  >  120*.  however,  the  picture  changes:  the  circumferential 
velocity  component  w^/u.  calculated  from  the  surface  pressures  is  higher  than  that  given 
by  inviscid-f low  theory  and  the  shape  of  the  resulting  distribution  is  considerably  dif¬ 
ferent  from  chat  given  either  by  inviscid  flow  or  that  obtained  by  measurements  at  the 
boundary-layer  edge.  The  result  of  these  differences  is.  as  shown  in  more  detail  below, 
that  if  the  boundary-layer  calculations  are  performed  for  the  external  velocity 
components  deduced  from  the  surface  pressures,  they  are  in  good  agreement  with  those 
computed  by  the  inviscid  input  data  so  long  as  4  is  less  than  120*.  For  4  >  120*  the 
results  differ  and.  perhaps  fortuitously,  the  inviscid  data  give  closer  agreement  with 
experiment. 

The  second  test  case  corresponds  to  the  Lockheed/NASA  Ames  C-wing  for  which  there 
is  experimental  data  obtained  In  the  NASA  Ames  6*  x  6'  transonic/supersonic  wind  tunnel 
with  a  semi-span  model  for  a  Mach  number  range  of  M.,  >  0.25  to  0.95  There  are  two 
series  of  tests,  one  for  the  measurement  of  surface  pressures  with  transition  free  and 
the  other  for  boundary-layer  measurements  with  transition  fixed  at  x/c  »  0.045  on  both 
upper  and  lower  surfaces.  The  boundary  layers  were  measured  only  at  two  chordvise  loca¬ 
tions  corresponding  to  x/c  >  0.218  and  0.421  at  the  mid-semi-span  location. 

In  contrast  to  the  first  test  case,  which  has  been  examined  by  all  of  the  boundary- 
layer  calculation  procedures  of  Chapter  4  and  by  the  Navier-Stokes  method  of  Chapter  5. 
the  second  test  case  has  been  calculated  only  by  the  Navier-Stokes  method  described  in 
Chapter  5  and  by  the  interactive  boundary- layer  method  described  later  in  this  chapter. 
As  a  consequence,  the  problem  of  the  specification  of  the  external  boundary  condition 
described  above  in  connection  with  the  first  test  case  does  not  arise. 

The  next  two  sections  are  concerned  with  Integral  and  field  methods.  In  Section 
6.1.  the  emphasis  is  on  the  prolate  spheroid  at  a  >  10*  with  some  results  for  the  higher 
angle  a  -  30*.  Section  6.2  presents  results  Cor  the  two  test  cases  referred  to  above  and 
provides  a  brief  description  of  the  interactive  procedure. 


6.1  Integral  methods 
P.D.  Smith,  RAE,  Farnborough,  Hants,  UK 


6.1.1  The  ellipsoid  at  10°  incidence 

This  is  the  only  test  case  to  which  the  integral  methods  described  in  Chapter  6  have  been  applied. 
For  each  method  results  are  shown  far  the  standard  test  case  using  the  external  flow  data  given  by  Stock 
together  with  some  examples  of  the  attempts  which  have  been  made  to  improve  the  agreement  between  calcula¬ 
tion  and  experiment. 

For  the  RAE  method,  Smith  has  produced  results  both  for  the  basic  method  and  for  the  method  using 
Cross's  velocity  profiles.  The  use  of  the  more  sophisticated  velocity  profiles  produces  a  small  improve¬ 
ment  in  the  results  but  large  differences  between  calculation  and  experiment  remain  on  the  leeward  side 
(4  >  120®  say)  of  the  body. 
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For  the  ONERA/CERT  method,  Gleyzes  notes  that  the  external  flow  directions  are  given  by  Stock  do 
not  agree  with  the  measurements  particularly  over  the  rear  leeward  side  of  the  body.  Gleyzes  has,  there¬ 
fore,  recalculated  the  external  flow  by  matching  it  to  the  experimental  data  at  XQ/2a  ■  0.48  .  As  will  be 
seen  from  the  figures  this  improves  the  agreement  between  experiment  and  calculation  but  significant 
differences  remain. 

For  the  ONERA  method,  Le  Balleur  and  Lazareff  have  calculated  both  the  standard  test  case  and  one 
in  which  the  external  flow  is  assumed  to  be  the  inviscid  potential  flow  which  would  be  present  in  the 
absence  of  all  viscous  effects.  This  second  calculation  improves  the  agreement  with  experiment  and  under¬ 
lines  the  unexpected  sensitivity  of  the  predictions  to  the  distribution  of  the  external  flow  angle.  This 
type  of  calculation  has  also  been  made  by  both  Smith  and  Gleyzes. 

As  will  be  seen  from  the  figures  all  of  the  methods  produce  results  of  a  similar  quality  and  for  the 
standard  test  case  they  agree  well  with  experiment  over  the  windward  side  of  the  body  but  on  the  leeward 
side  they  agree  more  closely  with  each  other  than  they  do  with  experiment. 
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6.2  Field  methods 


T.  Cebecif  Douglas  Aircraft,  CA,  USA 


6.2.1  The  ellipsoid  at  10°  incidence 

The  calculations  for  this  test  case  were  performed  by  Cebeci  tDACJ.  Hirschel  and 
Monnoyer  [MBB],  Humphreys  (FAA],  Menter  [DPVLRJ.  and  van  den  Berg  et  al.  [NLRJ  using  the 
boundary- layer  procedures  of  Chapter  4b  with  external  velocity  distributions  which  will 
be  indicated. 

The  results  obtained  by  Cebeci  are  described  in  detail  in  [3]  and  made  use  of  the 
external  boundary  condition  determined  from  the  solution  of  inviscid  flow  equations.  The 
distributions  of  momentum  and  displacement  thickness  show  similar  trends  to  those  of 
Section  6.1  where,  in  general,  the  external  boundary  conditions  determined  by  Stock  were 
used,  but  here  particularly  for  ^  greater  than  120*  the  calculations  are  in  much  better 
agreement  with  measurements.  There  is  a  tendency,  for  angles  of  120®.  for  the  calcula¬ 
tions  to  underpredict  at  low  values  of  x/2a  and  to  overpredict  at  the  furthest  downstream 
station.  This  is  probably  due  to  the  external  velocity  boundary  condition,  as  described 
in  [3]  and  discussed  above  in  relation  to  Fig.  1.  The  skin-friction  results  display  the 
experimental  trends  and  are  within  16%  of  the  measured  values  with  some  underprediction 
in  regions  which  are  consistent  with  the  experimental  pressure  distribution  as  argued  in 
[3].  The  flow  angles  show  even  better  agreement  except,  again,  in  the  downstream  region 
where  the  boundary  conditions  ace  subject  to  uncertainty. 

The  results  obtained  by  Hirschel  and  Monnoyer  [MBB]  made  use  of  the  inviscid  external 
velocity  distribution  and  a  brief  description  of  the  method  and  the  results  provided  by 
Monnoyer  ace  given  below. 

The  calculations  were  performed  on  the  1:6  ellipsoid  at  10®  incidence.  The  Reynolds 
number  based  on  the  ellipsoid  major  axis  was  7.700.000  and  the  freestceam  Mach  number 
0.17.  Transition  was  tripped  at  x/2a  »  0.2. 

The  mesh  had  89  nodes  in  the  axis  direction,  61  points  equally  spaced  In  the  cross- 
flow  direction  (half  body),  and  61  across  the  boundary  layer.  The  Initial  profiles  were 
computed  at  x/2a  »  0.0013,  and  the  calculation  was  laminar  until  x/2a  >  o.2:  intermlt- 
tency  was  applied  until  x/2a  *  0.226,  where  the  flow  became  fully  turbulent.  The  outer 
flowfield  was  obtained  from  the  potential  flow  solution,  and  no  displacement  effect  was 
applied  to  the  inviscid  pact  of  the  flow. 

Results  of  the  calculation  neglecting  (MBB-1)  and  including  (MBB-2)  the  effect  of 
surface  curvature  are  compared  with  experiment.  The  circumferential  distributions  of  the 
pseudo  two-dimensional  displacement  thicknesses,  the  flow  angles  at  the  wall  and  the 
skin-friction  coefficients  ace  plotted  at  four  stations.  In  addition,  the  wall  pressure 
distribution  is  provided.  The  comparison  of  both  calculations  shows  that  curvature  plays 
a  role  on  the  boundary-layer  flow  as  the  boundary-layer  thickness  increases,  namely  on 
the  leeside  region.  Since  the  velocity  profiles  at  the  vicinity  of  the  wall  are  only 
marginally  modified,  no  significant  changes  ace  to  be  seen  on  the  skin  friction  and  the 
wall  flow  direction,  although  there  is  a  trend  towards  the  measured  values.  On  the  other 
hand,  the  effect  of  curvature  on  the  metric  coefficients  away  from  the  wall  largely 
affects  the  solution  at  the  outer  pact  of  the  boundary  layer  as  illustrated  by  the  dis¬ 
placement  thicknesses. 

Concerning  the  comparison  with  experiment,  only  wall  flow  angles  Yw  agree  satisfac¬ 
torily.  except  on  the  leeside  at  the  last  station  where  it  is  clear  that  the  analytical 
potential  flow  solution  does  not  hold  anymore.  It  is  interesting  to  note  that  much  bet¬ 
ter  agreement  between  the  calculated  and  measured  values  is  obtained  if  the  latter  are 
slightly  "translated"  as  follows:  the  flow  angles  at  each  station  are  shifted  so  that 
the  symmetry  condition  (Yw  *  0)  at  the  leeside  is  satisfied  leading  to  shiftings  not 
larger  than  0.7®;  the  skin-friction  coefficients  are  translated  downwards  by  a  constant 
value  of  Acf  «  0.00025,  in  order  to  fit  them  with  the  computed  ones;  a  shift  Acp  •  0.025 
is  applied  to  the  measured  pressure  so  that  better  agreement  with  the  potential  fiow 
solution  is  obtained,  especially  on  the  windward  side.  it  is  indeed  surprising  that 
discrepancies  ace  smaller  at  the  leeside  where  it  is  expected  that  the  deviation  of  the 
actual  pressure  distribution  from  the  theoretical  inviscid  pressure  should  be  larger  than 
at  the  windward  side. 

The  results  obtained  by  Humphreys  made  use  of  the  experimental  pressure  distribution 
as  extracted  by  Stock.  They  show  Increasing  deviations  from  the  measurements  as  the 
circumferential  angle  increases  beyond  120*  with  some  discontinuities  which  may  stem  from 
the  details  of  the  numerical  procedure. 

Menter  [DFVLR]  has  used  a  second-order  accurate  finite-difference  method  which  is 
based  on  the  standard  and  the  zig-zag  box  method  (4-5].  He  has  performed  four  sets  of 
computations  by  using  different  boundary  conditions  at  the  boundary-layer  edge.  A  des¬ 
cription  of  the  results  provided  by  him  is  given  below. 
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For  Che  first  sec.  Che  boundary  conditions  were  taken  fron  potential  flow  theory. 
Figure  al  shows  the  wall  streanllne  angles  and  Fig.  a2  the  local  displaceaent  thickness 
at  four  different  x-stations.  The  conputations  ace  generally  in  good  agreenent  with  the 
expeciBents  except  in  a  small  region  in  the  rear  part  of  the  body  at  ^  z  140*.  There  a 
discontinuous  (weak)  solution  builds  up  in  the  direct  mode  computations.  This  discon¬ 
tinuity  is  certainly  a  feature  of  the  boundary-layer  eguations  itself  [6]  when  solved  in 
the  direct  mode  and  could  not  be  removed  by  changing  the  grid  resolution. 

A  second  set  of  calculations  has  been  performed  by  using  the  external  velocity  field 
computed  by  Haase  and  Stock  [7]  on  the  basis  of  the  experimental  pressure  distribution. 
An  improvement  of  the  results  could  not  be  achieved,  as  can  be  seen  in  Figs,  bl  -  b3. 
The  reasons  are  first  that  the  computations  are  still  performed  in  the  direct  mode  so 
that  singularities  in  the  solution  are  not  excluded  and  second  the  computations  of  [7] 
do  not  automatically  fulfill  the  symmetry  condition  at  the  leeward  line  of  symmetry, 
leading  to  the  unphysical  behavior  of  the  boundary- layer  method. 

The  singular  behavior  can  be  avoided  by  using  inverse  boundary  conditions  as  shown 
by  [0-9].  In  the  present  case  half  inverse  boundary  conditions  were  used  by  specifying 
first  the  length  of  the  velocity  vector  at  the  boundary-layer  edge  from  the  experimental 
wall  pressure  distribution  via  the  Bernoulli  eguation  and  second  the  direction  of  the 
wall  streamline.  This  experimental  input  was  used,  because  it  can  be  measured  quite 
accurately.  The  direction  of  the  external  streamline  was  computed  with  a  Newton  itera¬ 
tion.  Figures  cl  -  c3  show  the  results  (remember  that  cl  shows  the  input).  One  can  see 
that  the  discontinuity  in  the  solution  is  avoided  and  that  the  results  (Fig.  a2)  are  in 
good  agreement  with  the  measurements. 


In  order  to  avoid  the  need  of  an  experimental  input,  first  results  were  produced  by 
coupling  the  Euler  and  the  boundary- layer  equations.  The  boundary  layer  was  simulated 
with  the  wall  transpiration  approach  [10]  in  the  space-marching  parabolized  Euler  solver. 
The  correction  for  the  wall  pressure  was  computed  at  the  boundary-layer  edge,  defined  as 
follows: 


at  y  -  & 


Figures  dl  -  d3  show  the  results  of  this  viscous-inviscid  interaction  method.  The 
singularity  in  the  solution  is  completely  avoided,  as  was  expected  and  the  computations 
are  In  much  better  agreement  with  the  measurements  than  those  of  Figs,  al  -  a3.  It 
should  be  noted  that  the  major  shortcoming  of  the  viscous-inviscid  Interaction  approach 
is  the  arbitrariness  in  the  definition  of  the  boundary- layer  edge  (coupling  location). 
This  problem  can  only  be  avoided  by  using  parabolized  or  full  Naviei-Stokes  methods 
[11,12]. 

Calculations  performed  by  J.  I.  van  den  Berg.  Lindhout  and  B.  van  den  Berg  of  NLB. 
for  the  prolate  spheroid  at  10*  angle  of  attack  were  carried  out  with  BOLA.  the  NLR 
calculation  system  for  three-dimensional  boundary  layers.  The  calculations  were  started 
at  Xo/2a  «  0.30  with  given  initial  viscous  flow  conditions.  As  a  boundary  condition 
the  surface  pressure  distribution  was  employed,  the  flow  angle  at  the  boundary- layer 
edge  being  obtained  from  an  integration  of  the  Euler  equations.  A  description  of  their 
results  provided  by  them  is  given  below. 

In  a  first  report  [13]  grid  dependence  and  the  effect  of  changing  the  empirical 
constants  in  the  turbulence  model  were  investigated.  Refinement  of  the  grid  was  found 
to  have  little  influence  and  the  choice  of  slightly  different  constants  in  the  turbulence 
model  did  not  lead  to  essentially  different  results,  in  a  second  report  [14],  the  influ¬ 
ence  of  initial  data  and  boundary  conditions  on  the  calculation  results  was  investigated. 
Three  calculations  were  performed  subject  to  the  following  Initial  and  boundary 
conditions: 

1.  Measured  pressure  distribution:  Initial  flow  conditions  at  Xo/2a  *  0.30  from  Cebeci 
(based  on  boundary-layer  calculations  from  the  nose  using  the  potential  flow  bound¬ 
ary  condition) . 

2.  Potential  flow  pressure  distribution:  initial  flow  conditions  at  Xo/2a  •  0.30  from 
Cebeci . 

3.  Measured  pressure  distributions:  initial  flow  conditions  at  Xo/2a  *  0.30  from 
Gleyzes  (different  edge  stream  angles,  obtained  by  fitting  calculations  to  available 
measurement  data  at  x^/Za  -  0.40). 

The  first  case  represents  the  initial  and  boundary  conditions  as  defined  for  this 
test  case  by  the  Working  Group. 

The  calculation  results  for  the  three  cases  and  the  experimental  data  at  the  sec¬ 
tions  Xo/2a  •  0.48.  0.56.  0.64  and  0.73  are  plotted  in  a  series  of  graphs  appended. 
For  the  test  case  conditions  defined  by  the  Working  Group  (case  1).  the  calculation 
results  are  seen  to  compare  not  well  with  experiment  at  the  leeward  side  of  the  prolate 
spheroid.  The  region  of  flow  convergence  is  predicted  essentially  too  close  to  the  sym¬ 
metry  line.  The  calculations  with  the  same  initial  conditions  and  the  potential  flow 
pressure  distribution  (case  2)  produce  results  in  much  better  agreement  with  experiment, 
though  the  displacement  thickness  peak  is  substantially  overestimated.  Comparatively 
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SBooth  results  are  obtained*  probably  because  the  initial  and  boundary  conditions  used 
are  snooth  and  compatible.  Best  overall  agreement  is  obtained  when  the  measured  pressure 
distribution  is  used  together  with  initial  edge  flow  angles  derived  indirectly  from  the 
measurements  (case  3). 

It  is  concluded  that  the  test  case  is  dominated  to  a  great  extent  by  the  boundary 
and  initial  conditions.  More  specifically  boundary  and  initial  conditions  can  probably 
not  be  prescribed  independently.  The  poor  results  for  case  1  seem  to  be  due  to 
prescribing  the  measured  pressure  distribution  together  with  initial  edge  stream  angles 
from  potential  flow.  The  solution  of  the  Euler  equations  downstream  of  the  initial  line 
will  lead  then  to  an  essentially  wrong  edge  stream  angle  distribution. 

The  objective  of  the  Working  Group  is  to  determine  the  limits  of  boundary- layer 
calculation  procedures  in  3-D  separation  regions.  The  last  figure  shows  the  computed 
boundary-layer  thickness  distribution  at  two  sections  in  the  rear  part  of  the  prolate 
spheroid.  No  problems  were  encountered  in  performing  the  boundary- layer  calculations, 
but  the  large  circumferential  thickness  variations  indicate  that  the  boundary- layer 
approximation  (comparatively  large  derivatives  normal  to  the  surface)  does  not  hold 
anymore  and  determines  the  limit  here. 

6.2.2  The  C-wing 

The  calculations  were  performed  at  DhC  by  Kaups,  Mehta  and  Cebeci  [15]  at  the  design 
angle  of  attack  of  S  degrees,  at  Mach  numbers  of  0.70.  0.82,  0.8S  and  0.90  for  a  Reynolds 
number  of  6.8  x  10®.  The  interactive  boundary-layer  approach,  based  on  the  strip-theory 
approach  to  three-dimensional  flow,  has  been  described  in  [16]  and  incorporates  viscous 
effects  into  the  inviscid  flow  through  a  surface-blowing  boundary  condition.  Two 
inviscid-f low  procedures  corresponding  to  a  full-potential  code  [17]  and  an  Euler  code 
were  used  in  the  calculations.  The  transonic  potential  code  was  used  in  nonconservative 
form  with  a  numerical  grid  of  161  x  25  x  33  mesh  points  and  the  Euler  code  145  x  25  x  31 
mesh  points  arranged  in  a  C-grid  in  the  streamwise  and  spanwise  directions,  in  some  cases 
with  an  H-grid  in  the  spanwise  direction. 

Figure  2  allows  comparison  of  the  measured  and  calculated  results  at  -  0.70.  The 
interactive  calculations  with  the  Euler  code  and  H-grid  spanwise.  and  with  full  poten¬ 
tial  code  agree  well  with  experimental  data  except  near  the  leading  edge  where  the  flow 
is  supercritical.  As  expected,  the  calculated  lift  coefficients  also  agree  well  with 
the  normal  force  coefficient  -  0.483  calculated  from  the  measured  pressure 
distributions.  It  Is  evident  that  calculations  with  the  Euler  code  and  C-grid  in  both 
directions  reproduce  the  velocity  peaks  more  accurately  except  near  the  wing  tip  where 
the  local  lift  is  overpredicted  over  the  whole  chord.  The  results  of  Figure  3  correspond 
to  «  0.82  and  show  attached  flow  with  shocks.  The  interactive  calculations  with  the 
Euler  code  are  closer  to  the  measurements  than  those  with  the  full-potential  flow  solu¬ 
tions.  The  C-giid  was  used  in  both  directions  mainly  to  improve  results  near  the  wing 
tip  but  it  is  clear  that  better  agreement  with  experiments  has  resulted  near  the  shock 
location  and  over  the  whole  wing  apart  from  the  wing  tip  where  the  lift  tends  to  be  over¬ 
predicted,  and  caused  the  wing  lift  coefficient  to  exceed  the  measured  normal  force  coef¬ 
ficient  s  0.S3.  It  should  be  noted  that  the  predicted  pressure  recovery  near  the 
trailing  edge  at  the  wing  tip  differs  in  all  cases  from  the  measured  values  Indicating  a 
strong  decambecing  effect  in  the  experiment. 

Figure  4  shows  results  at  the  design  condition  at  which  flow  visualization  has  shown 
the  flow  in  the  outer  30  percent  of  the  span  to  be  separated  due  to  strong  shock- 
boundary-layer  interaction.  The  calculations  are  in  good  agreement  with  experiment 
except  in  the  region  of  the  separation  bubble*  and  the  corresponding  lift  coefficient  is 
only  slightly  lower  than  the  measured  value  of  of  0.54.  It  is  surprising  that  the 
results  obtained  by  interaction  with  the  full  potential  flow  method  are  in  better  agree¬ 
ment  with  those  with  the  Euler  method  than  the  results  of  Figure  3.  The  pressure  distri¬ 
butions  obtained  at  Ma  «  0.9  and  the  full-potential  code  are  again  in  reasonable  agree¬ 
ment  with  experiments  on  the  inboard  and  middle  portions  of  the  wing,  as  shown  In  Figure 
5:  the  predictions  on  the  outboard  portion  of  the  wing  show  attached  flow  whereas  the 
experiments  suggest  separation  because  the  predicted  shock  pressure  rise  is  spread  out 
over  a  considerable  distance.  Despite  this  discrepancy,  the  predicted  lift  coefficient 
is  close  to  the  measured  normal  force  coefficient  of  0.56. 
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Figure  1.  Comparison  of  calculated,  measured,  and  extracted  streamwise  and  circumferential 
velocity  components  for  the  proUted  spheroid  at  (a)  x/2a  *  O.Att,  (b)  x/2a  *  0.73 
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Fig.  a3:  Skin  friction  coefficient  at  four  different  x-stations  in 

comparison  with  experiments  (potential  pressure  distribution) 
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fa3:  Local  (fisplacement  thicknesses  at  four  different  x-stations 
in  comparison  with  experiments  (boundary  layer  edge  velocity 
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Fig.  c3:  Skin  friction  coefficient  at  four  different  x-stations  in 
comparison  with  experiments  (Inverse  boundary  conditions) 
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.  d2:  Local  displacement  thicknesses  at  four  different  x-stations 
in  comparison  with  experiments  (with  viscous-inviscid  inter¬ 
action) 
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Figure  3.  Comparison  of  calculated,and  experimental  chordwise  pressure  distributions  for 
=  0.82,  Re  =  6.8  x  10°.  (a)  Euler  with  H-grid.  (b)  Euler  with  C-grid. 


CHAPTER  7.  CONCLUSIONS 


The  present  study  of  the  Working  Group  10  was  confined  to  turbulent  flows  using 
both  the  field  and  integral  boundary  layer  equations.  Two  test  cases  were  selected  for 
numerical  simulation,  i.e.  the  DFVLR  prolate  spheroid  and  the  NASA-Ames  Wing  C. 

In  order  to  give  a  brief  summary  of  the  results  achieved  some  concluding  remarks 
are  listed  below  as  stated  by  the  various  WG  members: 

-  For  the  10*^  prolate  spheroid  case,  the  boundary  layer  equations  are  valid,  except  in 
the  regions  of  strong  convergence  of  the  wall  streamlines. 

The  10°  case  is  very  sensitive  to  azimuthal  distribution  of  the  given  external  flow. 

-  Problems  of  definition  of  external  flow  lead  to  lack  of  agreement  with  experiment  at  10°. 
Results  very  much  dependent  of  edge  conditions  even  at  10°. 

Boundary  layer  results  for  30°  agree  better  with  Navier-Stokes  results,  perhaps  because 
edge  conditions  and  integral  quantities  are  coherent  Iwhich  is  not  evident  for  10^). 

-  Wing  C  complicated  by  interference  from  root. 

3-D  boundary  layer  assumptions  are  very  sensitive  to  external  flow  angles. 

-  Apply  boundary  layer  theory  only  where  it  is  valid,  strong  and  global  interaction  should 
be  handled  by  coupling  approaches. 

For  attached  flow  or  for  flow  with  mild  cross  flow  separation,  the  viscous-inviscid 
interaction  method  is  accurate  and  efficient  and  well  suited  for  the  design  process. 

This  is  true  for  integral  as  well  as  for  field  methods. 

It  seems  that  boundary  layer  calculations  can  go  fairly  far  if  proper  boundary  con¬ 
ditions  are  introduced. 

The  initial  aim  was  too  ambitious. 

From  the  last  remark  one  can  easy  recommend  a  follow-up  Working  Group  which  would 
be  useful  if  viscous/inviscid  interaction  methods  are  available.  Furthermore,  it  is  re¬ 
commended  to  work  on  basic  feature  of  turbulent  flow  such  as  transition  and  turbulence 
modelling. 
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Appendix 

NLR  CALCULATIONS 

J.l.  van  den  Berg,  J.P.F.Lmdhout  and  B.  van  den  Berg,  NationaJ  Aerospace  Laboratory. 
NLR,  Anthony  Fokkerweg  2. 1059  CM  Amsterdam.  The  Netherlands 


Calculations  have  been  performed  for  one  test  case:  The  prolate  spheroid  at  10^  angle  of  attack.  The 

calculations  have  been  carried  out  with  BOLA,  the  NLR  calculation  system  for  three-dimensional  boundary 

layers.  The  calculations  were  started  at  x  /2a  ■  0.30  with  given  initial  viscous  flow  conditions.  As  a 
boundary  condition  the  surface  pressure  dis’iribution  was  employed,  the  flow  angle  at  the  boundary  layer 
edge  being  obtained  from  an  integration  of  the  Euler  equations. 

In  a  first  report  {J.l.  van  den  Berg  et  al,  19d7)  grid  dependence  and  the  effect  of  changing  the 

empirical  constants  in  the  turbulence  model  were  Investigated.  Refinement  of  the  grid  was  found  to  have 
little  influence  and  the  choice  of  slightly  different  constants  in  the  turbulence  model  did  not  lead  to 
essentially  different  results.  In  a  second  report  (J.l.  van  den  Berg  et  al,  1988)  the  influence  of  initial 
data  and  boundary  conditions  on  the  calculation  results  was  investigated.  Three  calculations  were  performed 
and  the  results  of  these  will  be  given  in  this  appendix.  The  Initial  and  boundary  conditions  for  the  three 
cases  are: 

1)  Measured  pressure  distribution;  initial  flow  conditions  at  x  /2a  •  0.30  from  Cebecl  (based  on 
boundary  layer  calculations  from  the  nose  using  the  potential  fl(&  boundary  condition) . 

2)  Potential  flow  pressure  distribution;  initial  flow  conditions  at  x  /2a  •  0.30  from  Cebecl. 

0 

3)  Measured  pressure  distribution;  initial  flow  conditions  at  x  /2a  ■0.30  from  Gleyzes  (different  edge 

stream  angles,  obtained  by  fitting  calculations  to  available  Wasurement  data  at  x  /2a  •  0.48). 

0 

The  first  case  represents  the  initial  and  boundary  conditions  as  defined  for  this  test  case  by  the  Working 
Group. 


The  calculation  results  for  the  three  cases  and  the  experimental  data  at  the  sections  x  /2a  •  0.48, 
0.56,  0.64  and  0.73  are  plotted  In  a  series  of  graphs  appended.  For  the  test  case  conditions  defined  by  the 
Working  Group  (case  1),  the  calculation  results  are  seen  to  compare  not  well  with  experiment  at  the  leeward 
side  of  the  prolate  spheroid.  The  region  of  flow  convergence  is  predicted  essentially  too  close  to  the 
symmetry  line.  The  calculations  with  the  same  Initial  conditions  and  the  potential  flow  pressure 
distribution  (case  2)  produce  results  In  much  better  agreement  with  experiment,  though  the  displacement 
thickness  peak  Is  substantially  overestimated.  Cotaparatlvely  smooth  results  axe  obtained,  probably  because 
the  initial  and  boundary  conditions  used  are  smooch  and  compatible.  Best  overall  agreement  is  obtained  when 
the  measured  pressure  distribution  is  used  together  with  Initial  edge  flow  angles  derived  Indirectly  from 
Che  measurements  (case  3). 

It  Is  concluded  chat  the  test  case  is  dominated  to  a  great  extent  by  Che  boundary  and  Initial 
conditions.  More  specifically  boundary  and  initial  conditions  can  probably  not  be  prescribed  Independently. 
The  poor  results  for  case  1  seem  to  be  due  to  prescribing  the  measured  pressure  distribution  together  with 
initial  edge  stream  angles  from  potential  flow.  The  solution  of  the  Euler  equations  downstream  of  the 
initial  line  will  lead  then  to  an  essentially  wrong  edge  stream  angle  distribution. 

The  objective  of  the  Working  Group  Is  to  determine  the  limits  of  boundary  layer  calculation 
procedures  in  3D  separation  regions.  The  last  figure  shows  the  computed  boundary  layer  thickness 
distribution  St  two  sections  in  the  rear  part  of  the  prolate  spheroid.  No  problems  were  encountered  in 
performing  the  boundary  layer  calculations,  but  the  large  circumferential  thickness  variations  indicate 
that  the  boundary  layer  approximation  (comparatively  large  derivatives  normal  to  the  surface)  does  not  hold 
anymore  and  determines  the  Limit  here. 
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